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Abstract 

This  work  is  based,  in  part,  upon  the  Kubo  transform 
technique  for  calculating  transport  coefficients  by  way  of 
space-time  correlatJon  functions.   These  space-time 
correlation  functions  are  evaluated  by  a  special  path  integral 
formulation  first  described  in  a  paper  by  Percus  and  Yevick. 
This  integral  is  the  basis  of  a  minimum  principle  which 
results  in  bounds  upon  space-time  correlations.   The  path 
integral  representation  is  used  to  derive  various  auto- 
correlation functions,  e.g.  entropy-entropy;  moreover,  a 
method  for  representing  these  auto-correlations,  via  a 
minimum  prin.ciple  is  presented.   In  order  to  establish 
a  usable  minimum  principle  the  integrand  ha?  to  be  guaranteed 
positive  definite.   This  is  accomplished  by  an  extension 
into  the  imaginary  time  domain.   Thus,  with  analyticity  in 
the  upper  half  complex  time  plane  and  sufficient  fall  off 
at  infinity  a  lower  bound  inequality  for  transport  coefficients 
can  be  given. 

As  a  simple  test  case  the  idea  is  applied  to  a  model 
chain  of  oscillators  and  a  second  such  system  with  different 
parameters.   A  simple  exact  inequality  for  the  self-diffusion 
coefficient  is  obtained.   A  direct  comparison  between  an 
ordinary  perturbation  expansion,  to  first   order,  for  a 
potential  of  the  quartic  interaction  type  and  a  perturbation. 


to  like  order,  in  the  path  integral  formulation,  for  the 
same  interaction,  is  made.   The  two  methods  are  seen  to 
give  the  same  result  following  a  renormalization  procedure. 
That  is,  the  presence  of  certain  oscillating  terms  in  the 
integrand  of  the  path  integral  formulation  are  eliminated 
by  renormalizing,  about  a  frequency  shifted  harmonic 
oscillator  term. 

Finally  the  path  integral  representation  is  used  to 
discuss  the  N-body  problem  of  transport  in  a  Debye-type  crystal 
with  nearest  neighbor  cubic  forces.   To  guarantee  normaliza- 
tion, one  starts  with  a  harmonic  system,  for  which  an  exact 
solution  is  known  and  turns  on  the  remaining  quartic 
potential  slowly  from  time  minus  infinity  to  time  zero  via 
a  decaying  exponential.   The  density  subject  to  Liouville's 
equation  has  its  normalization  miaintained  since  the  partial 
time  derivative  of  the  phase-space  integral  of  the  density 
vanishes.   One  is  then  able  at  time  equal  to  t  to  inquire 
as  to  the  correlation  of  the  one-body  phase-space  density 
with  that  at  time  equal  to  zero,  aliC'Wing  the  calculation 
of  transport  coefficients.   Specifically  an  inequality  for 
the  velocity-velocity  correlation  function  is  determined  and 
its  Infinite  imaginary  time  limit  is  examined. 
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CHAPTER  I 
Theory  of  Transport  Coefficients  and  Time 
Correlation  Functions:   A  Survey 

Pick's  law  of  diffusion  J   =  -  DVn ,  Ohm's  law  of 

P 

electrical  condution  J   =  -aVA   and  Fourier's  law  of 

e 

thermal  conduction  q  =  -  AVT   express  the  generalized  forces 
required  to  maintain  a  small  steady  flux  of  particles, 
charge,  and  heat  respectively.   They  may  also  be  used  to 
describe  the  approach  to  equilibrium  of  a  system  initially 
in  disequilibrium  due  to  the  presence  of  gradients  in 
concentration  n ,  electrostatic  potential  4),  or  temperature  T. 
The  constants  of  proportionality  are  termed  linear  transport 
coefficients:   diffusion  D,  electrical  conductivity  a, 
thermal  conductivity  X.      We  shall  be  concerned  with  the 
derivation  of  these  transport  coefficients,  and  their 
generalizations,  from  knowledge  of  the  microscopic  properties 
of  the  system  involved.   Let  us  first  survey  the  approaches 
to  thio  problem  which  have  been  used  heretofore. 

A.   Boltzmann  Equation  Method. 

The  Boltzmann  equation  describes  the  nonequilibrium 
behavior  of  the  single-particle  distribution  function  of  a 
dilute  monatomic  gas  subject  only  to  binary  collisions. 
In  order  to  obtain  the  transport  coefficients  one  assumes 
that  the  syster,  is  close  to  local  equilibrium;  this  is  a 
state  for  which  equilibrium  is  established  in  volumes  which 
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contain  large  numbers  of  particles  but  are  still  small 
on  the  scale  of  macroscopic  variations .   In  each  volume 
region  one  has  local  values  of  temperature  T,  density  ri,  and 
mean  velocity  u.   By  expressing  the  single-particle 
distribution  function  as  a  local  equilibrium  distribution 
plus  a  correc'-ion  proportional  to  the  gradients  in  the  local 
equilibrium  quantities  one  can  obtain  a  solvable  linearized 
version  of  the  Boltzmann  equation.   The  method  of  solution, 
after  Chapman  and  Enskog,  depends  upon  the  assumption  that 
the  single-particle  distribution  function  is  time  dependent 
only  through  the  time  variation  of  the  local  equilibrium 
variables.   Thus,  one  obtains  the  normal  solutions  presumed 
valid  for  a  system  close  to  equilibrium,  and  also  explicit 
expressions  for  the  transport  coefficients.   Prom  the 
Boltzmann  equation  it  is  possible  to  generalize  the  macro- 
scopic laws,  e.g.  Fourier's  law,  by  considering  corrections 
to  the  linearized  Boltzmann  equation  and  including  nonlinear 
terms.   Typically,  Fourier's  law  becomes  q  =  -  XVT  + 
+(9V*uVT  +  other  terms  of  like  order)  +  higher  order  terms. 
The  method  of  H.  Grad  is  to  expand  the  distribution 
function  in  terms  of  Hermite  polynomials,  the  coefficients 
of  which  are  taken  as  state  variables.   The  number  of 
coefficients  taken  depends  on  the  degree  of  approximation 
desired,  and  the  complete  set  of  Hermite  coefficients  is 
equivalent  to  the  distribution  function.   If  the  expansion 
is  stopped  after  a  finite  number  of  terms,  the  distribution 


function  Is  obtained  as  a  function  of  a  number  of  Hermite 
coefficients;  consequently,  a  determined  system  of  differ- 
ential equations  can  be  obtained  for  the  coefficients. 
By  selecting  a  sufficient  number  of  Hermite  coefficients 
rather  general  solutions  of  the  Boltzmannequation  can  be 
approximated. 

The  first  systematic  method  for  obtaining  all  the 
higher  order  density  corrections  to  the  Boltzmann  equation 
was  accomplished  by  Bogoliubov.   Basic  to  his  approach  was 
that  each  higher  density  correction  to  the  Boltzmann  equation 
takes  into  account  the  interaction  of  more  particles. 
However,  Bogoliubov's  method  depends  upon  the  functional 
assumption  which  states  that  the  distribution  functions  for 
two  or  more  particles  do  not  depend  explicitly  upon  the 
time.   Instead,  these  functions  depend  upon  time  only 
by  way  of  the  single-particle  distribution  function.   More 
recently  M.  S.  Green  gave  an  alternative  derivation  of  the 
generalized  Boltzmann  equation  without  recourse  to  the 
functional  assumption.   Green's  work  is  based  upon  cluster 
expansion  methods. 

B.   The  Liouville  Equation  Approach. 

The  method  of  time  correlation  functions  is  an 
attempt  to  base  a  linear  macroscopic  transport  equation  theory 
directly  on  the  Liouville  equation.   The  theory  owes  its 
origin  to  Einsten  and  Nyquist  and  of  late  to  M.  S.  Green, 
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Kubo,  and  Mori.   The  essential  Idea  of  the  time  correlation 
function  method  Is  to  consider  the  effect  of  a  local 
equilibrium  state  on  the  complete  N-partlcle  distribution 
function  for  the  system.   All  of  the  obtainable  Information 
about  the  system  Is  contained  In  the  N-partlcle  distribution 
function.   In  the  time  correlation  function  method  It  Is 
assurred  that  the  N-partlcle  distribution  function  can  be 
written  as  a  local  equilibrium  N-partlcle  distribution 
function  plus  correction  terms.   The  local  equilibrium 
function  depends  upon  the  local  macroscopic  variables, 
tem^perature ,  density  and  mean  velocity  and  upon  the  position 
and  momenta  of  the  N  particles  in  the  system.   The  correction 
to  this  distribution  function  is  determined  via  Liouvllle's 
equation  and  the  assumption  that  at  some  initial  time  the 
system  was  in  local  equilibrium  but  at  a  later  time  is 
tending  towards   complete  equilibrium.   Hence,  as  expected, 
the  normal  solutions  to  Liouvllle's  equation  can  be  found 
and  an  exp-ression  for  the  corrections  to  local  equilibrium 
in  powers  of  the  gradients  of  the  local  variables  as  well. 
By  keeping  only  the  first  term  in  the  gradient  expansion 
one  derives  the  generalized  linear  macroscopic  transport 
equations.   These  expressions  depend  upon  the  dynamics 
of  all  N  particles  in  the  system  and  apply  to  any  system, 
regardless  of  its  density. 
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C.   The  Kubo  Transform  and  Related  Material. 

The  formulas  for  transport  coefficients  usually  have 
the  form  of  a  one-sided  Fourier  transform  of  a  time  correla- 
tion function  (defined  as  the  relevant  ensemble  average  of 
a  two-time  product  A(t)B(t')). 


(1.1)  a(uj)  = 


e~^"^<A(0)B(t)>  dt  . 


0 

The  variable  w  in  the  transform  is  to  be  regarded  as  a 
frequency  and  if  necessary  to  provide  convergence  it  can 
have  a  small  imaginary  part.   In  most  cases  only  the  low 
frequency  behavior  is  of  interest.   At  first  thought  it 
appears  that  in  order  to  calculate  a  ti-ansport  coefficient 
by  means  of  time-correlations,  one  has  to  carry  out  the 
horrid  task  of  solving  the  general  N-body  problem.  However, 
it  is  really  not  that  bad;  for  example,  one  needs  only  the 
time  dependence  of  certain  quantities  A(p,q)  and  this 
information  often  is  more  accessible  than  the  complete 
solution  of  the  equations  of  motion. 

The  quantum  mechanical  expressions  for  transport 
coefficients  are  more  complicated  than  their  corresponding 
classical  counterparts.   As  R.  J.  Zwanzlg  points  out  the 
following  changes  are  to  be  made.  (1)  The  functions  A  are 
replaced  by  their  corresponding  quantum  mechanical  operators 
(2)  The  equilibrium  phase  space  distribution  function  is 
replaced  by  the  equilibrium  density  matrix  p.  (3)  The 
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ensemble  average  is  calculated  by  taking  the  trace  of  the 
operator  product,  <A>  =  Trace  (pA).   (4)  The  time  derivative 
of  the  operator  A  is  given  by  Heisenberg's  equation  of 
motion   dA/dt  =  i[H,A]  where  H  is  the  system's  Hamiltonian 
operator.   The  solution  of  this  equation  of  motion  is 
A(t)  =  e^^"(A)  e~^^^.   (5)   Lastly,  the  defining  relation 
for  the  time  correlation  function  is  altered  by  the  use 
of  the  Kubo  transform.   The  Kubo  transform  of  an  operator 
A  is  another  operator  A  defined  by 

6 


(1.2) 


~         1    ,,   XH/.N   —  XH 
A  =  g-   dX  e   (A)  e 

0 


in  which  B  =  1/kT,  where  k  is  Boltzmann's  constant  and 

T  is  the  absolute  temperature.   Thus  transport  coefficients 

in  quantum  mechanics  have  the  form 


(1.3) 


a(co)  = 


e"^'^^<  A(0)B(t)>  dt 


0 


At  sufficiently  high  temperatures  the  Kubo  transform  of  an 
operator  is  just  the  operator  itself  and  thus  the  quantum 
mechanical  expressions  reduce  back  to  their  classical 
counterparts . 

D.   Kubo's  Method 

Theoretical  connections  between  time-correlation 
functions  and  transport  coefficients  have  been  studied  from 
several  points  of  view.   A  distinction  is  usually  made 


-6- 


between  transport  processes  driven  by  external  forces  and 
those  associated  with  Internal  forces.   The  response  of  an 
equilibrium  system  to  an  externally  Imposed  force  lends 
Itself  to  a  methiod  after  Kubo.   Kubo's  method  Is  simple  and 
operational.   For  Instance,  consider  the  Kubo  method  applied 
to  electrical  conductivity.    Basically  a  measurement  of 
electrical  conductivity  consists  In  the  observation  of  the 
electric  current  responding  to  an  applied  electric  field. 
Inasmuch  as  the  experimental  measurement  Is  carried  out 
repeatedly,  one  must  average  the  results.   The  conductivity 
Is  then  defined  by 


(1.4) 


<r>  =  aE  +  O(E^) 


where  <J>  Is  the  average  electlrc  current,  a  the  conductivity 
and  E  the  applied  field. 

To  measure  the  frequency  dependence  of  the  electrical 
conductivity  one  turns  on  a  time-dependent  electric  field 
E(t)  at  t  =  0  and  then  measures  the  average  current  <J(t)>. 
One  finds 


(1.5) 


<J(t)>  = 


dt'  $(t-t' )E(t' )  +  O(E^) 


0 


where  $(t-t')  is  an  after-effect  function  showing  the 
response  at  t  to  a  disturbance  at  t ' .   One  defines  by 
Fourier  analysis  frequency  components  of  the  current  and 
field. 
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(1.6) 


<J^>  =  I  e"^^^  <J(t)>  dt 
0 

oo 

E    =  f  e"^'^^   E(t)   dt 
-to     J         ~ 

0 


By  rearrangement  one  has 


(1.7) 

where 


<J  >  =  a(a))  E 
~(i)      '  '  -co 


a(oo)  = 


^  e-^'^^  Ht)    dt  . 


0 


Kubo's  statistical  mechanical  theory  is  merely  a  model 
of  this  experiment.   As  a  useful  aside  one  may  derive, 
following  Kubo,  the  classical  response  of  a  system  due  to 
a  dynamical  disturbance.   The  system  Hamiltonian  is  assumed 
to  be  of  the  form 


H  =  H„  +  H   , 
0    ext 


H   ,  =  -  A  F(t) 
ext 


where  F(t)  is  a  force  and  H„  is  the  unperturbed  Hamiltonian. 
The  statistical  distribution  function  of  the  system  under 
discussion  follows  the  equation  of  motion 


9t 


=    -I    ( 


3f^  aH, 


^^t  ^"t. 


9q   9p     8p   3q 


=  -^\'h^ 


If  the  external  disturbance  is  small  the  distribution 
function  may  be  expressed  as  an  unperturbed  part  plus  a 
small  perturbation  term. 


-8- 


f^  =  f  +  Af 


where  f  satisfies 


9f 
9t 


=  (H,f)  =  0  . 


By  algebraic  manipulation  one  determines  Af 


9Af 
dt 


=    (H,Af)  -  F(t)(A,f) 


This  can  be  rewritten  as 


9Af 

at 


=  iLAf  -  F(t) (A,f ) 


in  which  iLg  =  (H,g)  is  the  Polsson  bracket  operator, 
The  natural  motion  of  any  dynamical  quantity  A(t)  is 
given  by 


A(t)  =  A(P^)  -  e"^""^  A^^t-T^ 


where  P  is  a  phase  point  and  exp  (-ixL)  represents  the 
natural  unperturbed  motion  of  the  system.   This  equation 
for  A(t)  is  evidently  a  solution,  corresponding  to  the 
Heiser.berg  picture  in  quantum  mechanics,  of 


dA 
dt 


=  (A,H)  . 


The  solution  of  the  Poisson  bracket  operator  equation  is 

t 

Af(t)  =     gi(t-t')L  ^^^^^    ^^,.  F(t') 

—  oo 

since  F(-°°)  =  0,   AF(-°°)  =  0;   that  is,  the  system  is  in 
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equilibrium  at  t  =  -<».   The  expectation  value  for  any 
quantity  B  is 


<AB(t)>  = 


B  Af(t)  dr 


dt 


dr  B  e^^^-^'^^  (A,f)  F(t') 


_00 

t 


dt'  F(t' ) 


(A,f )  B(t-t' )  dr 


for  which   dr   is  a  volume  element  in  phase-space. 
Thus,  rewriting 

t 

<AB(t)>  =   J  <D^^(t-t')  F(t')  dt' 

—  00 

with  the  response  (after-effect)  function  given  by 


$g^(t)  =  I  (f,A)  B(t)  dr 


f(A,  B(t))  dr 
If  the  distribution  function  f  is  canonical 


f  =  c  e 


-3H 


and  therefore. 


(f,A)  -  3Af 


which  finally  allows  one  to  write 
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-J-g^Ct)  =  6<A(0)B(t)> 


=  -  B<A(0)  B(t)> 

since  ^  <A(t),  B(t+T)>  =  0  . 

Returning  to  the  model  conductivity  experiment 
it  can,  in  like  manner,  be  shown  [17]  that 

(1.8)  "^ba^^^  "  ^^'^a'^b^^^^ 

in  which  $,   is  the  after-effect  function  giving,  in  this 
D  a 

case,  the  current  in  the  b-th  direction  due  to  a  field  in 
the  a-th  direction  and  thus  $(t)  is  a  time  correlation 
function.   The  corresponding  frequency  dependent  electric 
conductivity  tensor  is 

-itot 


(1.9)       ^ba^""^  "  ^ 


<J  J.(t)>  dt 
a  b 


0 


Transport  coefficients  associated  with  internal 
disturbances  e.g.  thermodynamics  can't  be  handled  directly 
by  Kubo's  formulation;  nevertheless,  other  approaches  do 
exist  and  among  them  are  the  indirect  Kubo  method,  the 
Fokker-Planck  technique,  the  regression  hypothesis  and  local 
equilibrium  theory.   Each  of  these  approaches  will  be 
briefly  reviewed. 
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E.   The  Indirect  Kubo  Method. 

In  order  to  clarify  the  indirect  Kubo  method  a 
hydrodynamic  problem  will  be  discussed.   At  t  =  0  an 
external  force   F(R,t)   is  applied  to  a  fluid  in  thermal 
equilibrium.   This  force  causes  fluid  motion  and  gives  rise 
to  a  velocity  field   v(R,t).   Expressing  the  relation  between 
force  and  response  as  Fourier  components  one  has 


(1.10) 


F(aj,k)  = 


v((jo,k)  = 


F(R,t)  e-i^-?-'^'^^  dt  d^R 


,-o   <_  \   -ik'R+itot  ,,  ,3t, 
v(R,t)  e   ~  ~      dt  d  R 


f 
The  solution   of  the  hydrodynamic  equations  of  motion  if 

the  force  is  transverse  i.e.  V*F  =  0  and  VxF  -^   0  leads  to 


(1.11) 


v(co,k)  =  a(co,k)  F(a),k) 


a  transverse  velocity  field.   Hence  the  coefficient  a(co,k) 

is 


tt 


(1.12) 


2        -1 
a(co,k)  =  [nk   -  icop]    , 


tt 


for  forces  small  enough  so  that  the  equations  of  motion 
may  be  linearized  in  deviations  from  equilibrium. 

the  calculation  which  leads  to  this  equation  is  macroscopic 
and  follows  entirely  from  the  Navier-Stokes  equation 
which  is  discussed  briefly  in  Chapter  VII. 
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where  p  is  the  mass  density,  and  n  the  shear  viscosity 
coefficient.   To  obtain  the  shear  viscosity  coefficient 
from  a(aj,k)  one  goes  to  the  limit  of  long  wavelength  viz. 
k  ^  0  and  then  to  the  limit  of  low  frequency: 

2  2 

(1.13)  n  =  lim  lim   ^-^  Re  a(a3,k)  . 

03^0  k->0   k 

However,  a(a),k)  can  be  obtained  from  Kubo's  method 
since  the  force  P(R,t)  is  externally  applied  and  v(R,t) 
is  a  response  to  that  force.   Hence,  Kubo's  method  leads  to 
a  statistical  mechanical  formula  for  a(a),k)  and  consequently 
n  as  well. 

P.   The  Fokker-Planck  Technique. 

The  Pokker-Planck  theory  is  based  upon  ideas  like 
those  occurring  in  Brownian  motion,  assuming  that  the  motion 
is  determined  by  a  large  number  of  small  impulses.   Let 
P(Rj^,  tg/R,  t)  be  the  probability  per  unit  volume  of  finding 
the  particle  at  position  R  at  time  t,  given  its  position  Rq 
at  initial  time  t^.   This  function  PCRq'  ^n^^'  ^^  ^^  determined 
by  way  of  the  Fokker-Planck  equation  [10] 

(1.14)  W  ^   ^'  ^-^^  "^  V^(DP)  . 

To  Interpret  A  and  D  let  AR  be  the  displacement  of  the 
particle  during  a  time  interval  At  which  is  small  in  a 
macroscopic  sense  but  larger  than  characteric  molecular  times; 


+ 

the  duration  of  a  molecular  collision  is  of  the  order  of 

-12 
10     seconds. 
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then  one  finds  in  deriving  the  Fokker-Planck  equation  that 
A  is  an  average  velocity. 


(1.15) 
and  that 
(1.16) 


A  =  lim 
At->0 


D  =   lim 
At^O 


<AR> 
At 


<(AR)^> 
^At 


The  formula  given  for  D,  the  coefficient  of  self- 
diffusion,  was  originally  deduced  by  Einstein.   Hence,  we 
have  in  principle  a  methxod  of  obtaining  a  transport 
coefficient  by  the  use  of  a  time  correlation  function. 
If  one  writes  the  displacement  of  the  particle  as  the  time 
integral  of  the  velocity 


At 


(1.17) 

then  D  is  given  by 


AR  = 


v(£)  an 


At  At 


(1.18)   D  =   lim 
At^O 


1 
FAt 


<v(2.)  •v(5,'  )>  d£  d£' 


0   0 

where  the  integrand  contains  the  auto-correlation  function 
of  the  particle  velocity.   By  translation  in  time  invariance 
one  may  write 


(1.19) 


<v(i)-v(£')>  =  <v(0) •v(£'-£)> 


and  thus  integrating  one  gets 


(1.20)    D  =   1 


im 


At-^0 


At 
1    ■ 

3   J 
0 


(1-  ^)  <v(0)-v(£)>  d£ 
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The  term  l/^t   may  be  left  out  provided  At  is  small  but 
nonetheless  macroscopic,  and  if  the  auto-correlation 
function  decays  fast  enough: 

At 


(1.21) 


D.  1 


<v(0) •v(£)>  d£ 


0 


G.   The  Regression  Hypothesis. 

This  method  is  based  upon  Onsager's  hypothesis 
concerning  the  regression  of  fluctuations.   The  hypothesis 
states  that  spontaneous  fluctuations  from  thermal  equilibrium 
decay,  on  the  average,  according  to  the  macroscopic  transport 
laws.   Kubo,  Yokota,  and  Nakajima  proposed  a  method  based 
upon  this  hypothesis.   This  is  illustrated  as  follows. 
Consider  the  macroscopic  transport  law  concerning  the 
average  value  a(t)  of  a  dynamical  quantity  a(t). 


(1.22) 


da(t)  ^ 
dt 


M  a(t) 


where  clearly  M  is  a  rate  constant  i.e.  transport  coefficient 
for  the  decay  of  a(t).   Multiply  both  sides  of  the  transport 
law  by  the  initial  value  a(0),  and  then  average  over  an 
equilibrium  ensemble  of  initial  values  obtaining 


(1.23) 


<^   5(t)  a(0)>  = 


-  M<a(t)  a(0)> 


If  aCt)  is  the  numerical  value  of  a  function  A(p(t),q(t)) 
in  phase-space,  then  one  may  write  the  average  in  the  form 
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of  a  time-correlation  function. 


(1.24) 


<a(t)a(0)>  -  <A(p(t),q(t))  A(p(0),q(0))> 


=  <A(t)A(0)> 


The  regression  hypothesis  has  been  used  implicitly. 
That  is,  in  averaging  over  an  equilibrium  distribution  of 
initial  a(0),  the  major  contribution  comes  from  small  values 
of  a(0),  on  the  order  of  spontaneous  fluctuations,  but  one 
wants  the  transport  laws  to  be  valid  for  just  such  small 
initial  values . 

Following  the  above  equations  one  may  write  the 
transport  law  as 


(1.25) 


<A(t)A(0)>  =  ~  lvi<A(t)A(0)> 


Integrating   this  equation  over  the  small  time  interval  (0,t) 

t 


(1.26) 


t 


0 


<A(t ' )A(0)>  dt '  =  -  M 


<A(t' )A(0)>  dt 


0 


and  if  Mt  <<  1  then  one  obtains 


(1.27) 


<A(t ' )A(0)>dt '    -  -Mt<A(0)A(0)>+  O(M^t^) 


0 


Using  time  invariance  to  shift  the  time  origin  i.e. 


(1.28) 


<A(t '  )A(0)>  =  <A(0)A(-t ' )> 


and  introducting  another  time  derivative 
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(1.29) 


A(-t')  =  A(0)  - 


A(-Jl)  di 


^ives  for  t  small, 

t 

<A(0)A(0)>dt '  - 


f 


t  t  • 


(1.30) 


0 


<A(0)A(-£)>d£  dt' 


0  0 
=  -  Mt<A(0)A(0)>   . 


1  dA 
The  first  term  on  the  left  vanishes  since  <AA>  =  <2   gT~-* 

viz,  the  equilibrium  average  of  a  rate  of  change  must 

vanish.   In  the  second  term  the  time  origin  is  shifted 

back  by  an  interval  I   thus  replacing  <A(0)A(-£,)>  by 

<A(il)A(0)>  and  lastly  the  order  of  integration  is  reversed, 

Hence,  performing  the  t'  integration  and  dividing  both 

sides  of  the  equation  by  t  gives, 

t 
(1.31)       I  (1-  ^)  <A(£)A(0)>  dt  =  M<A^>  . 

0 

Provided  that  the  time  correlation  decays  fast  enough  one 
may  in  analogy  with  the  Fokker-Planck  technique  drop  the 
£/t  factor  and  obtain 

T 


(1.32) 


M  = 


<A(£)A(0)>  /  <A  >  d£ 


0 


where  the  upper  limit  has  been  replaced  by  a  time  t  lon^ 
enough  for  the  integrand  to  have  decayed  to  zero. 
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H.   The  Local  Equilibrium  Theory 


This  method  was  devised  by  Mori  and  the  basis  of  it 
is  the  observation  that  when  a  typical  nonequllibrium  system 
relaxes  toward  thermal  equilibrium  it  passes  through  a  series 
of  local  equilibrium  states.   Consider  the  case  where  only 
a  single  macroscopic  variable  A  is  of  interest.   The  ensemble 
mean  value  of  this  quantity  is  denoted  by  A.   In  an  equilibrium 
ensemble  where  A  is  held  fixed  by  some  kind  of  external 
constraint,  the  phase-space  distribution  function  f(p,q)  is 


(1.33) 


f(p,q)  = 


^-SH-yA 


where  the  partition  function  0(3, y) 


(1.34) 


:(3,y)  - 


^-BH-yA 


dp  dq 


is  the  normalizing  factor  and  H  is  the  system  Hamiltonian. 

This  distribution  function  gives  maximum  entropy  under 
the  constraints  of  fixed  H  and  A.   The  normal  thermal 
equilibrium  distribution  function,  appropriate  when  A  is 
not  constrained,  is  obtained  by  setting  y  =  0-   The 
I'elationship  between  the  Lagrange  multipliers  3  and  y ,    and 
the  mean  values  H  and  A  is  given  by 


(1.35) 


H  =  - 


9  lo£ 


93 


A 


9  lof 


9y 


The  main  assumption  of  the  local  equilibrium  theory 
is  that  the  phase-space  distribution  function  maintains  the 
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above  form  for  all  times.   However,  the  assumption  is  not 
correct,  except  in  the  trivial  case  where  A  like  H  is  a 
constant  of  motion.   Moreover,  it  has  the  additional  failing 
that  the  mean  value  of  the  rate  of  change  of  A  vanishes. 

Nevertheless,  the  true  distribution  function  of  a 
relaxing  system  probably  doesn't  deviate  appreciably  from 
the  locally  canonical  form  and  therefore  the  ideas  of  local 
equilibrium  theory  are  used  in  practice.   One  allows  a  system 
to  grow  old  for  a  certain  time  interval  t  which  is  long 
compared  to  molecular  collision  times  but  short  enough  such 
that  the  macroscopic  state  of  the  system  hasn't  altered 
drastically.   Then  the  mean  value  of  the  curret  J  =  A 
is  measured  at  x.   The  average  is 


(1.36)      J(T)  =  Q  ^ 


!  ,    V   -3H-yA  ,   -, 
A(t)  e     '   dp  dq 


Small  deviations  from  equilibrium  correspond  to  small  Yj 
and  the  thermodynamic  force  associated  with  A  is 

(1.37)  F  =  ky  +  O(y^) 

where  k  is  the  usual  Boltzmann  constant.   One  then  expands 
the  average  current  in  powers  of  y 

(1.38)  J(t)=  <J(t)>-y<A(t)A(0)>  +  y<J(t)><A(0)>  +  O(y^) 

and  since  <J(t)>  is  the  current  at  thermal  equilibrium  it 
vanishes,  so  that 

(1.39)  J(t)  =  -  y<A(t)A(0)>  +  O(y^)  . 
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Following  the  manipulative  procedure  of  the  already  outlined 
regression  hypothesis,  equation  (1.39)  becomes 


(1.40) 


J(t)  =  L  F(0) 


where  the  transport  coefficient  L  is  given  by 

T 


(1.41) 


L  = 


k 


f 


0 


<A(£)A(0)>  dil 


The  various  methods  reviewed  all  lead  to  the  same 
results  with  the  following  reservation.   Some  methods  lead 
immediately  to  expressions  for  transport  coefficients  involving 
integrations  over  characteristic  time  intervals  x,  long 
compared  with  molecular  times,  and  short  compared  with 
macroscopic  times .   Other  methods  lead  to  integrals  over  an 
infinite  time  interval  together  with  some  limiting  process 
(e.g.  (jo  ->-  0 )  .   For  all  practical  purposes,  even  though 
different  formally,  these  two  basic  approaches  should  be 
considered  identical. 
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CHAPTER  II 

A.   Introduction  and  Basic  Formulation 

In  equilibrium  statistical  mechanics  all  of  the 
thermodynamic  averages  of  a  system  are  known  if  its 
partition  function  in  an  external  field  can  be  calculated. 
Similarly  for  the  space-time  correlations  required  in  the 
theory  of  transport  processes,  it  is  possible  to  Introduce 
a  special  time-dependent  field  and  a  time-dependent  partition 
function  to  achieve  the  same  effect.   By  adding  a  time- 
dependent  inhomogeneous  term  to  be  described  one  may  define 
the  time-dependent  partition  function  as  the  complete 
phase-space  integral  of  a  suitable  N-body  many  time  Boltzmann 


factor,  Z  = 


t 


G(x,p)   1  fdxCt')  dp(t')  (where  x(t')  and 
~  ~   t'=0   ~      ~ 
p(t')  are  the  possible  positions  and  momenta  for  any  time 

interval  between  0  and  t).   To  obtain  the  weight  function 

for  expectations  which  involve  any  x(t)  and  p(t)  (0  <_  t  <_  t )  , 

one  must  guarantee  that  x(t)  and  pd)  be  restricted  by  the 

equations  of  motion  and  Initial  conditions.   Now  the  equations 

of  motion  can  be  guaranteed  by  applying  the  appropriate  delta 

functions;  whereupon, 

(2.1)  G(x,p|0,t)  =  TT  6(x(t)-  ^  p(t))6(p(t)+  i|ill)G[x  (0  )  ,p  (0  )  ] 

T=0 

for  which  G[x(0),p(0)]  is  the  initial  unnormallzed  distribu- 
tion usually  to  be  taken  as  an  equilibrium  Boltzmann  factor. 
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After  Fourier  representation  of  the  delta  functions,  one 
obtains  as  the  required  generall2'.atlon  (to  within  normalization) 
the  path  Integral  [l8] 


G(x,p|0,t)  = 
t 


exp  [1 


(x(t)-  i  p(t))-^(t)  dx] 


(2.2)     -exp  [1  I  (p(t)+  i|ill).n(T)  dT] 

0 

•  D{?(t)}  D{n(T)}  G[x(0),p(0)]   . 

This  is  to  be  the  basis  of  a  variational  principle  and 
in  fact  a  minimum  principle  resulting  in  bounds  upon  space- 
time  correlations.   As  prototype,  we  consider  the  case  of 
a  simple  harmonic  oscillator.   The  path  Integral  representa- 
tion will  be  used  to  derive  the  one-body,  one-dlmenslonal 
velocity  and  entropy  auto-correlation  functions.   In  addition, 
a  method  for  normalizing  this  simple  problem  will  be 
presented . 

B.   Velocity  Auto-Correlation  for  Simple  Harmonic  Oscillator 


Before  considering  the  path  Integral  for  this 
one-body  example  a  correctly  normalized  expression  for  the 
velocity  auto-correlation  will  be  derived  by  using  a  simple 
phase-space  Integral.   The  harmonic  oscillator  solution  in 
terms  of  initial  values  is 


x(t)  =  x(0)  cos  tjo„t  + 

0     CO 


x(0) 


sin  a)„t 
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where  the  phase-space  density  is 
p(x(0),x(0))  =  -^^   e 


3niWQ   -  -^(x^CO  +  oaQX^CO)) 


subject  to  I  p(x(0),x(0))  dx(0)  dx(0)  =  1. 


Therefore  one  may  write 


R 


-   ^  x(t)x(0) 
<e  > 


3maj, 
~2Tr 


I-       r      -   f  R(x    (O)cos    0Jnt  + 
e 


^0' 


(2.3) 


_oo     _oo 


,       X(0)X(0)  .  .x  r>  •o  OO 

+    -S;;^ ^"^    '^O^)       -    ^(x2(0)+a.2x2 


(0)) 


dx(0)    dx(0) 


After  integrating  and  differentiating  with  respect  to  R 
at  R  =  0  one  obtains 


(2.4) 


<x(t)x(0)> 


Bmco 


p—  cos  (J0„t  . 


For  time  translational  invariant  systems 


<x(t)x(0)>  =    <x(t+t' )x(t' )> 


and  therefore 


=  -  <x(t+t ' )x(t ' )>  by  parts  integration 


=  -  <x(t)x(0)>      as  t'  ^  0 


3t' 


<x(t)x(0)> 


(2.5) 


<x(t)x(0)>  =  o—  COS  u^t 

3m      0 
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Returning  to  the  unnormallzed  path  integral 
representation  one  has  in  one-dimensional  form 


oo        oo 


<G(x(t),x(0),x(t),x(0))>= 


G(x(t),x(0),x(t),x(0)) 


_oo        —CO 


t        .      .        t ' 
-  imBx^(O)   -6(j)(x(0))   i  /  ^  (t '  )  [x  (t  '  )-x  (0  )-  i  /  F  (x(t  "  )  )dt"  ]dt ' 

e  e  e  '^  o 


t  t- 

fd^Ct')  ]  rdx(t')  dx(0)  dx(t)  dx(0)  , 


0 


0  + 


where  the  initial  velocity  distribution  has  been  taken  as 
Gaussian  and 


G  =  e 


-  i  Rx(t)x(0) 


with 


(})(x(0))  =  I  ma3QX^(0) 


F(x(t"))  = 


-  mco„x(t") 


then  once  again  the  problem  is  to  determine 

1 


<x(t)x(0)>  =  -  2  1^  <e   2 


t 


R  x(t)x(0) 


R=0 


By  integrating  out  the  x(0)    dependence  we  get 


See  Appendix  B, 
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Rx(t)x(0) 


1   ^22 


<e 


>  = 


m 


Ba)Qx"(0) 


■exp  - 


2m3 


(2.6) 


5(t'  )C(t")  dt'  dt"-  iRx(t) 


i  /  ^(t')x(t')dt' 


1   n2  '  2  /  ,  ^  n  O 

jj-  R  X  (t)]  •  e 


?(t')  dt' 


t  t  ' 
1  /  /  WQ^Ct' )x(t")  dt'  dt" 

e   °  °  D[ 


.] 


For  a  Gaussian  integral,  the  steepest  descent  method 
Is  identical  with  the  exact  result.   Hence  the  only  problem 
is  to  find  the  saddle  "point"  (function)  of  the  exponent. 
If  the  kernel  of  (2.6)  is  varied  with  respect  to  ^(x)  and  x(t) 
(with  T  =  t',t")  one  obtains  the  integro-dif ferential 
equations 


(2.7)   - 


m6 


C(t')  dt'  +  if^^  +  ix(T)  +  i 


WQx(t")  dt"  =  0 


(2.8) 


-  i  ?(t)  +  i  1  Wg  C(t')  dt'  =  0 


T-dif ferentlating  these  equations  and  re-substituting  the 
results  gives 
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?(t)  =  0 


a  boundary  condition,  and 


2 ' 
5(0)  =  i3m(jo„x  (0 )   which  is  a  coupling  condition 

between  equation  sets. 


The  results  of  the  x-dif f erentlation  are 


x(t)-0JqX(t)  =  0    subject  to  x(0)  =  x(0);  x(t)  =  x(t) 


and 


C(t)  -  Wq  5(t)  =  0 


Therefore,  solving  these  equations  subject  to  all 
four  conditions  yields 


(2.9)    x(t)  = 


x(t)  sin  oj^T  +  x(0)  sin  u)p(t-T) 
sin  (jo„t 


(2.10) 


2   x(t)-x(0)  cos  oa^t 

5(t)  =  16ma3Q( 2 )  ^°^    ^^(t-T) 

sin  co„t 


iRa)„   x(t)  cos  tO(-,t  -  x(0) 

+  — 2~(  2 ^  ^°^  tOQ(t-T) 

sin  oa^t 


Furthermore  the  kernel  can  be  simplified  if  equation  (2.7) 
is  multiplied  by  5(t)  and  integrated.   In  like  fashion  (2.8) 


is  multiplied  by  x(t)  and  also  integrated.   There  results 

t  t 

~  ^  J  J  ?(^)S(t')  dt'  dT  +  i  J  C(t)x(t)  dx 

0  0 

(2.11) 


IR 
2m3 


C(T)x(t)  dT  +  i 


t  T 


0  0 


cOq  C(T)x(t")  dt"  dT  =  0 
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t  t 


(2.12)   -  i 


5(t)x(t)  dT  +  i 


f   ? 

Wq  ?(t' )x(t)  dt'  dT  =  0 


0 


0  T 


Integrating  the  single  integral  in  (2.12)  by  parts  and 
adding  this  result  to  (2.11)  implies 

t  t 


C(t)x(t)  -  C(0)x(0)  =  - 


1_ 
m3 


5(T)C(t')  df  dT 


0  0 


(2.13)   + 


iR 
2m3 


^(T)x(t)  dT  +  i 


t  T 


0  0 


U)q  ^(T)x(t")  dt"  dT 


t  t 


+  i 


(      ? 

oJq  C(t'  )x(t)  dt'  dT  +  2i 


5(t)x(t)  dT 


0  T 
Thus  the  kernel  is 

1   o  2  2 


iR 


2 

1  ,   R   '2 


-|m6a.^x^(0)   ^x(t)  /  C(t')dt'   -^(-^x^(t)) 


•  e 


I  i(C(t)x(t)-C(0)x(0)) 


Therefore  equation  (2.6)  becomes 


-  I  Rx(t)x(0) 
<e  > 


oo      oo 

r 


_00     _00 


t 


-  |mBoaQX^(0)  ^x(t)/  ?(t')dt' 
e  e       o 


2 

^x^it)      i  i(at)x(t)-5(0)x(0)) 
e^""^       e^  dx(0)  dx(t)  , 


which  with  the  explicit  reintroduction  of  equations  (2.9) 
and  (2.10)  and  some  simple  algebraic  manipulations  becomes 
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(2.1^) 


CO     CO 


_oo  _00 


exp  -  i-(c-^^x^(0)  +  2c^2^^°)^^t) 


+  C22X  (t))  dx(0)  dx(t) 


where 


^11    ^22 


2    2 
3mw„+R(jj„  cos  oj^t 

-^   - 
sm  Wj-,t 


2  2        2 

3ma)|-,  cos  w^t    Rtjo„(l+  cos  w^t ) 


'12 


•  2   , 
sm  a)„t 


2  sm  cjo„t 


Simultaneously  integrating  this  quadratic  square 
and  differentiating  with  respect  to  R  evaluated  at  R  =  0 
yields  to  within  normalization 


27T       2 

(2.15)   <x(t)x(0)>  =  p— p-  a)„  cos  co„t  sin  Wj-jt 

(mB^Q) 


C.   Normalization, 


The  normalization  for  expression  (2.6)  is 


N  = 


-  p-  m3x  (0)  -  p-  m3a)„x  (0) 
e  e 


—  00       _oo 

t 


t' 
r 


(2.16)  -exp  [i    C(t ' ) (x(t' )-x(0)  + 

0  0 

t  t 

.fTdCCt')  T~Tdx(t')  dx(0)  , 
0  0 


coQx(t")dt")dt'] 


which  after  the  variations  are  carried  through  becomes 
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N    = 


f      -    p-  mpcOpX    (0 ) 
e 

_oo      _oo 


(2.17) 


2  2  2 

-,         „    X    (t  )-2x(0)x(t  )cos    0)   t  X    (O)cos    w„t 

-    2m3coQL 2 2 -' 

e  sin   oj^t  sin   oj^t 


Integrating 


•^0 


0 


dx(0)  dx(t) 


(2.18) 


N  = 


27T 


3m(jO 


p-  sin  a)„t  . 


If  expression  (2.15)  Is  now  divided  by  (2.l8)  we  get 


1 
<x(t)x(0)>  =  o—  cos  oa^t  , 


0 


which  Is  clearly  Identical  to  (2.5).   Aside  from  the 
fact  that  this  method  gives  the  correct  result  Its  drawback 
is  that  the  variational  principle  must  be  applied  separately 
to  both  the  expectation  of  a  function  and  its  purported 
normalization.   However,  this  objection  will  be  overcome 
in  the  full  N-body  problem  by  an  initial  normalization 
technique . 

D.   Entropy  Auto-Correlation. 


The  entropy  in  a  local  equilibrium  ensemble  is 

defined  by 

d 


where 


S(t)  =^I  (H.-h.)x. 


H.  -    contribution  of  the  j-th  molecule  to  the  Hamlltonlan 

h.  =  enthalpy  per  molecule 

X.  =  position  of  j-th  molecule. 
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For  simplicity  In  the  calculation  we  consider  only  the 
kinetic  energy  contribution  to  the  entropy  and  so 


S(t)  =  1^  I  n:x^(t)x(t) 


"  dt 


J(t) 


Because  of  time  translation  Invariance  the  auto-correlation 
may  be  expressed  as 

'  2 

<S(0)S(t)>  = ^  <J(0)J(t)> 

3t 


dt 


^  <|  mx^(0)x(0)  I  mx^(t)x(t)> 


If  one  chooses 


Q  ^  m^  x^(0)  e"^  x(0)x^(t)x(t) 


then  as  expected 


2 


<S(0)S(t)>  =   ^l^<^   x2(0)  e-«  x(0)x-(t)x(t)^ 


9f 


R=0 


2    2 
\   <^  x^(0)x(0)x2(^)^^^)> 


3t 
Thus ,    one  has 

<G(x(t),x(0),x(t),i(0))>= 


OO         00 


f   2  .  „ 

m    2 


TT 


xMO) 


_oo      _00 


-Rx(0)x2(t)x(t)   -  7   "^^^^(0)   -  I  m3(j2^2(0) 
•e  e  e 


1  /  ?(t')[x(f)-x(0)+  /   a)^x(t")dt"]dt' 
'e   D  o    ^  D[ . . .]  , 
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which  upon  completing  the  square  In  x(0)  and  Integrating  Is 


oo       oo 


^m6    '-m3 


\  5    2      -Rx(0)x^(t)x(t) 

at')  dt']^)  ^  e 


_oo      _00 


1   o  2  2 


exp[-  2  mgco^x  (0)-  ^ 


(2.20) 


t 


m6 


+  i    ?(t ' )x(t' )dt'  +  i 

0 
t  t- 


C(t' )^(t")dt'  dt" 


0 

t  t' 
f 


a3Q5(t')x(t")dt'  dt"] 


0  0 


~f  d?(t')l  fdxCt')  dx(0)  dx(t)  . 
0  0+ 


Using  the  same  techniques  as  In  expressions  (2.6)  to  (2.13) 
we  obtain 


<G(x(t), . . . )>  = 


OO      00 


_oo   _oo 


^m6    m3 


r 


0 


at')  dt']^) 


(2.21)   .  m!^-Rx(0)x2(t)x(t)  ^-|n^3.oV(0) 


i  l[at)x(t)-aO)x(0)] 
e  dx(0)  dx(t) 


where   x(t)   Is  given  by  (2.9)  and 


C(t)x(t)    -    ao)x(O)    =    IBitiWqC 


2  2  2 

P    X    (t  )-2x(0  )x(t  )cosaj„t  +  x    (O)cos    t)0„t 


sin   a3„t 


also 

4- 
f 


^  o    9    ii    x(t)-x(0)cos    co^t    5 

at')dt']^    =    B   m   (Oq( 2 ^^ 

sin'  Wj-jt 


)"[ 


cos    ajQ(t-t ')dt']' 


2  2  2 

P    P    P        X    (t  )-2x  (0  )x(t  )cosooQt+    X    (O)cos    to„t 

=    3   m  Wq    (    P ) 

sin   Wj^t 
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After  R-differentlating  (2.21)  evaluating  the  result  at 
R  =  0,  completing  the  square  In  x(t)  and  x(0),  integrating 
and  dividing  the  result  by  expression  (2.l8)  gives  the 
normalized 


9R 


<G(x(t),x(0),x(t),x(0))> 


R=0 


2    2  2  3 

m   a)„       -.         2    cos    w^t         cos    oo„t    -,  2oj„cos    co^t 

W~    ^m6  b\^  wl       ^mj  ~        sin   co^t 


^0 


^0 


3        3  3 

(jO^COS     Up,t  -.  2CO„COS     (JJp,t 

+   ^—. ^)^^  +   -^ —1    + 


(2.22) 


sin   w^t      2^2 
0      m   6 


•    2      , 
sm   w„t 


m^gS 


0 

2         3 
,         -4co„cos    a)„t 

^    [( ^ — 


sin    (ja„t 


2  2         5  2         2  2         5 

3(jO(-,    cos    Wj-jt         oj^cos    w^t         4a3„cos    cogt  2003^003    Wpt 


3 
sin   cjo„t 


sin   oj^t 


3 
sin   co„t 


sin   oa„t 


1^4  4        7  4        5 

4to„cos    Wp^t         2a)|-,cos    ajp.t         2a3p,cos    w^-jt 


5 
sin   (jjj-,t 


5 
sin   tOj-,t 


•    2      , 
;in   ojpjt 


To  obtain  <S(0)S(t)>  requires  taking  the  second 
time  derivative  of  (2.22)   but  as  this  is  of  no  Immediate 
consequence  it  has  not  been  done  but  only  mentioned  for 
completeness . 
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CHAPTER  III 
Imaginary  Time  Extension 

It  is  known  that  an  integral  of  an  exponential  always 
admits  to  a  variational  principle.   However,  in  the  present 
case  the  variational  principle  does  not  take  the  form  of 
the  more  useful  minimum  principle.   This  is  so  because  of 
the  oscillatory  nature  of  the  complex  exponential;  therefore, 
if  a  useful  principle  is  to  be  established  the  oscillation 
problem  must  be  eliminated.   This  will  be  done  by  extending 
all  the  time  factors  to  imaginary  times,  i.e.  t  ^  It  [18]. 
If  transforms  such  as  Cdx)  ^  15  (x)  are  employed,  then  it 
is  possible  to  guarantee  a  positive  integrand,  and  consequently 
a  rigorous  bounding  principle. 


We  will  require 


°? 


1 
2  J 


0 


f(t)  dt  for  even  functions  (i.e. 
f(t)  dt).   Assuming  analyticity  in  the  upper  half  complex 
time  plane  and  sufficient  fall  off  at  infinity,  we  have  the 
unique  imaginary  extension 


(3.1) 

It    follows    that 


f(iT)     =     ^ 


(    g^    g)    f(t)    dt    . 

t     +T 


f(t)    dt    =    lim   TTTfdx) 


lojt 


In  practice,  we  will  need  extension  to  [  f(t)  e"""""  dt 

—  00 

(if  only  to  allow  w  ->-  0+  and  therefore  eliminate  the  constant 
component).   Nevertheless  the  above  considerations  apply 
here  as  well: 
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(3.3) 


where 


f(t)  e^^^  dt  =  lim  -^xf  (ix 

T;->-oo 


(3.4) 


f  (It)  =  - 


:-5^)  e^""^  f(t)  dt 
t  +x 


A  concept  somewhat  similar  to  this  Is  used  In  the 
Quantum  statistical  mechanical  Green's  function  technique 
where  the  real  Green's  function  at  time  =  0  Is  matched 
with  Its  Imaginary  counterpart  at  time  =  -IB  (6  Is  reciprocal 
temperature  and  0  >  Im  (t-t*)  >  -3).   This  essentially 
matches  the  real  axis  with  a  line  through  the  Imaginary  axis, 
whereas  we  match  the  entire  real  line  with  a  contour  at 
Infinity  In  the  upper  half  plane. 

A.   Model  of  Time  Extenslor. 

As  a  demonstration  of  this  Idea  (for  w  ^  0+)   f(t)  Is 
taken  to  be  J|-.(t),  the  Bessel  function  of  the  first  kind. 


Thus, 

00 

(3.5)  f(t)  dt  = 

and  It  Is  known  [24]  that 
(3.6) 


jQ(t)  dt  =  2 


oo 

f 


2 — 2  Jq^^^-*  ^^   ^    [iQ(ax)  -  LQ(aT)] 


0 


t  +x 


with  a  >  0,   Re  X  >  0  and  Ip^(ax)  -   —  ]    cosh(ax  cos  6)  dG  Is 

o 
the  Bessel  function  of  the  second  kind,  while 

2   ^^2 
Lj-,(ax)  =  —    I   slnh(aT  cos  6)  dG  Is  the  modified  Sturve 


function.   Thus, 
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(3.7)   llm 


(3.8) 


IT 

r 


T  cosh(T  COS  0)de  -  2 


^/2 
r 


T  sinh(T  cos  9)  d9 


=  lim  2 


7T/2 


0 


T  e 


0 

-T  COS 


d0 


should  give  the  same  result  as  equation  (3.5).   But  expres- 
sion     (3.8)  can  also  be  written  as 


(3.9) 


tt/2 


lim 

2 

-[•->-oo 

lim 

2 

X->-oo 

T  e 


-T 


-T  sm 


de  = 


T  e     de  =   2  ; 


and  therefore  the  verification  is  made 
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CHAPTER  IV 

Introduction. 

We  can  now  set  up  the  required  variational  principle 

for  transport  coefficients,  each  of  which  can  be  written  in 
-B 


the  form 


X    e 


for  X  >  0  and  B  real.   This  is  done 


by  observing  that 


X(-B+l-ln  p)p  = 


6/6p 


,  -B    ,        -B 
Xe  at   p  =  e   , 

-B 


X(-B+l-  In  p)p  -  -  X(B+  In  p)  =  0  at  p  =  e 


and  6  /6p  6p 


X(-B+l-  In  p)p  =  -X/p  <  0  for  all  p  >  0 
-B 


Hence  setting  p  =  e   ,  we  have  the  following  lower  bound 
inequality 


(4.1) 


X  e 


-B 


X(l+B-B)  e  ^. . 


(one  obtains  equality  for  B  reducing  to  B).  This  enables  one 
to  make  comparisons  between  a  model  system  (B)  and  another 
system  (B)  gaining  information  about  the  latter.   As  a  simple 
test  case  this  idea  is  applied  to  a  model  one-dimensional 
chain  of  oscillators  and  a  second  such  system  with  different 
parameters.   A  simple  exact  inequality  for  the  diffusion 
coefficient  is  obtained. 

A.   Decoupling  the  Hamiltonian  [5,6] 


We  consider  the  Hamiltonian  (neglecting  end  effects) 

N  p         N 

H  =  I       (PV2m  +  cu  )  -  i  c  I       (u  u     +  u  u    ) 


j  =  l 

N 


J=l 


j=l 


>2 
J" 


=  I       (Pf/2m  +  cuf)  -    c      I      u,u 


2 

I . 

J 


N 


j=l 


J  J+1 
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The  first  term  describes  a  set  of  independent  harmonic 
oscillators,  one  on  each  site,  while  the  second  is  the 
interaction  between  each  oscillator  and  its  nearest  neighbor, 
By  introducing  the  canonical  transformation 


(4.3)    u,  =  4  I  (X 

k 


n    "   TtJ  L    v^i.  COS  kjd  -  r,  sin  kjd) 


'k 


(4.4)    p 


where 


j  =  71  A  ^"^"^k^k 


I    (moo  Xj^  sin  kjd  +  r,  cos  kjd) 


2c  1/2 

i^y.   =    [—  (1-  cos  kd)] 


d   =   inter-oscillator  distance 
subject  to  the  Poisson  bracket  rules 


we  get  after  length  manipulation 
(4.5) 


2      2 
,       r,    mu),   ^ 


k 


,t 


■2m 


2   "k' 


a  completely  decoupled   sum  of  ordinary  harmonic  oscillator 
Hamiltonians . 


B.   Simple  Inequalities 


The  diffusion  coefficient  can  be  written  via   t  ->-  It, 


in  the  form 


\    e 


-B 


and  therefore  the  inequality  holds 


the  usual  transformation  (see  Chapter  VI  and  Appendix  A) 
leads  to  a  Hamiltonaln  which  mixes  modes  of  wavenumber 
k  and  -k . 
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and  can  be  evaluated  in  any  convenient  form.   With  the 
Intention  of  computing  the  diffusion  coefficient  equation 
(4.1)  Is  made  explicit  by  selecting 


X  =  e 


-  i  I    R[x^(t)-x^(0)]^ 

k 


B  =  );  [|  mex^(0)+  I  Aj^Bx^(O)-  1 


(4.6) 


k 


-  x,  (0))dt'  -  - 
k  m 


?^(t')(x^(t') 


t  ' 
r 


C^(t')A^x^(t")dt'  dt"] 


0  0 


B   =    I    [^   m6x^(0)  +  J   A,  BxJ(O)-  1 


^k 


-  X,  (0))dt'  -  - 

k  m 


r 


^j^(t')(Xj^(t') 


0 


2  "k^"k 
t  t' 

C,  (t')A,  X,  (t")dt'  dt"]  , 


'k 


k  k 


0  0 


2      '     —2 
where  r,^(t)  =  mx,^(t),   A^  =  mco,  ,   A^  =  mco^ 


k 


k 


k 


As  previously,  we  seek  <[x,  (t  )-x,  (0 )  ]  >  - 

3     -  U  R[Xk(t)-x,._(0)]' 
=  -  2  ^^  <e    k 

^^  R=0 


so  as  to  ultimately 


establish  an  Inequality  for  the  diffusion  coefficient  based 
upon  Its  "model"  value  using  B.   We  will  act  as  If 
normalization  Is  automatic  —  It  Is  not  In  this  case  — 
and  show  later  how  it  can  be  made  so. 


(4.7) 


Thus,    the    right    hand   side    of    (4.1)    reduces    to 

t    t' 
'         2 


OO  00 

r 


_oo  _oo 


(1+    y    [4(A,  -A,  )Bx,   (0)+    - 
^2k      kk  m 


?^(f)(Aj^-A^)x^(t") 


0    0 


1      r.'2 


1    ^'2 


■dt'dt"])    exp    I    [-   |RxJ(t)+    Rx^(t)x^(0)+   |Rx^(0) 

is. 


\  mBxJ(O)-   I  k^   x^(0)+    1 


^^(t')(x^(t')-x^(0))    dt'    + 
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m 


t  t ' 
r 


0  0 


?j^(t' )A^x^(t")  dt'  dt"] 


t  t- 

.  TTdC(t')  "TTdx.Ct')  dx,(0)  dx  (t)  dx  (0)  . 

0  0+  K  K 

It  should  be  observed  that  this  kernel  Is  apart  from 
the  summation  over  index  k ,  essentially  that  dealt  with 
in  Chapter  II  and  so  by  applying  the  same  methods  we  can 
write  for  the  kernel  of  (4.7) 

(4.8)   exp  -  I  );  (c^^x^(0)+  2c^2^k^°^''k^^^'^  c^^x^it)) 

k 

where  the  x^(t)  are  given  by  (2.9),  the  ^^{t)    by  (2.10) 


and   the    c . .    are 

?         2 

Bma3^+Rco„cos    tij„t+   R 

c        =   — 2 -■ 

Rw„    cos    oj„t 


'22 


'k 


2         2 
3mco„+Roa„cos    a)„t    -R 

sin'''co„t 


?  2        2 

23ma)j-,cos  oo^t    RcOq(1+  cos  oj^jt) 

°i2  =     __...2.  ,       +        rr 


'0' 


sin  co„t 


sin  oj^t 


Hence,  expression  (4.7)   after  carrying  through  the 
lengthy  but  direct  algebraic  and  differential  operations  is 


Ir  I  A(B-B+l)e 


-B 


R-0    k 


=  I    ^a^   cos  03,  t- 


3m 


(4.9) 


1  ^  k  _  _j£-j  ^   /       k 


23m^m    m   /A     ^^ 

k 


)  , 


and  for  A  =  A^,  one  gets 


(4.10)     ^ 


X   e 


-B 


=  y  Q—  cos  oj,  t 


3m 


R=0    k 

Furthermore,  the  resultant  contribution  to  the  final 
inequality  is 
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3F  I  A  e 


R=0 


<    y  [75—  cos  03,  t 
—    ^    Km  It 


k 


3m 


(^.11) 


-.       A,    A,    / —   t  cos  OJ,  t 
1  /  k  _  Jkx  ■m   / k  N  -1 

2m 3  m    m  ^  J  A^  ^  sin  co^t  '  -^ 


where  the  sense  of  the  Inequality  has  changed  because  in 
taking  the  R-derivatlve  a  factor  of  -1  was  brought  across 
both  sides  of  the  equation.   By  combining  these  equations 
one  establishes  an  inequality  on  the  velocity  auto- 
correlations between  the  model  case  and  the  second  compara- 
tive system. 

2  M 


<[x^(t)-x^^(0)]">  -  <[xj^(t)-x^(0)]  > 


(4.12) 


1  I     23-A7  '^k^'^k-^k^^ -. —^ 

k     k  OJ,  sm  00,  t 


k     k 
and  perhaps  a  more  convenient  form  of  which  is 


<[x^(t)-Xj^(0)]2>  -  <[XjJt)-Xj^(0)]2>" 


(4.13) 


A 


k 


1  ^  2rA-  ^Af -1) 

k     k    k 


'       —2   2 
where    A,  =  A,  00,  /w. 


3      2 

CO,  t  cos  0),  t 

k 1^ 

sin  0),  t 


k    k  k'  k 


C.   Evaluation  of  Sums. 

In  order  to  evaluate  this  sum  a  transformation  to  an 
integral  form  will  be  made  where  for  a  system  of  N  one- 
dimensional  harmonic  oscillators  one  has  [19] 


-i|0- 


w. 


i  I  c< 


2 
N    ^  k  TT 


(•         cos    cot    dw 


/    2      2^1/2 

(cOj^-O)      ) 


for  ;;hlch  the  cutoff  frequency  Is 


mo)^  =  Ha        and   g(oj)  dco   = 


2  du) 
"T~^   2  s  1/2 


Is  the  fraction  of  circular  frequencies  between  oj  and 
0)  +  do),  and 


u 


L 


2 

TT 


cos    oot    do) 


,    2      2.1/2        "O'^^L 

0     ^'"l"'"  ^ 


=   J„  (w.  t )      as      N  ->■ 


Hence,  returning  to  Cl.lO),  we  can  write  with  the  introduc- 
tion of  the  cutoff  frequency  and  spectral  density 


1    '^L 

-   l,^   cos    (.j^t 


2  .  "L. 


2   '^L 


(^.14) 


COS  ojt 
/  2   2.1/2 

(cOj^-O)  ) 


do) 


Wt 


L 


Since  the  diffusion  coefficient  is  readily  obtained  by 
integrating  the  (absolutely  convergent)  velocity-velocity 
correlation,  we  find 


(^.15) 


D  =   lim 

03^0  + 


'  e^-^  <x(t)x(0)>  dt 


0 


~   2 

b 


CO 


L 


In  like  manner 


-Hi- 


N  f  23  A,   'A, 
k     k    k 


'        3     2 
A,       to),   cos  03,  t 


sin  oj,  t 
k 


(4.16) 


u. 


2 

7T 


'    (7.^-1) 


23  A,  'A, 
k   k 


3     2 
tw   cos  wt  dcjj 

~         ...  2   27172 
(sm  cot )  (oa^  -cjj  ) 


For  the  sake  of  simplification  In  Integration,  we  make 

the  assumption  that  A   and  A   are  Independent  of  site  and 

2       2 
with  the  use  of  the  trigonometric  Identity,  cos  cot  +sln  cot  =  1, 


consider  Instead 


0) 


L 


(4.17) 


'     -(^-1) 


23A  'A 


0 


2t       CO  /    ...     1   s  -I 

—  — ^ ■  (-  sm  cot+  — T— — 7r)doj. 

TT  /  2   2.1/2  sm  (jjt 

(co^-oa  ) 


In  order  to  establish  the  diffusion  Inequality,  we  break 
Integral  (4.17)  Into  two  parts  and  examine  first 


(4.l8)llm  ^i^  -1) 
a^O    ^^^  ^ 


OJ-r 


0    0 


2    or; 

Ti    f    2      2.1/2 


—at 
e  t    sin   ojt    dco   dt 


—at 
where  e     Is  the  usual  convergence  factor  [*  ]  and 


11m 
a^O 


e    t  sin  cot  dt 


0 


2^  2 
a  +C0 


^  0 


thus  (4.l8)  gives  a  zero  contribution  to  the  Inequality, 
At  this  point  let  the  Imaginary  time  extension  be  made 


t  ^  It  such  that  the  second  part  of  (4.17)   becomes 


COt 


(4.19) 


CO 


,    2      2sl/2      slnh    cot 


doj  dT 


See  J.  L.  Lebowltz  and  J.  K.  Percus,  The  Phys .  Rev.  155 
1  (1967). 
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After  carrying  through  the  Integration  of  ('4.19) 
one  gets 

(7 1)U   0)^  ) 


2eA  'A    "''2  "L' 

_2 

since  the  t  integral  gives  a  factor  of  u   and  this, 

3  2   2  -1/2 

in  turn,  reduces  the  u)   integral  to  co(cOy-a)  )     which 

yields  uij.      Thus,  we  have  established 
(^.20)  D  -  d'^  <  ^  (|^  -1) 


and  (4.20)  can  be  rewritten  (see  (4.15)) 


A        A'  -  A  ^A    ^ 


which  is  equivalent  to  the  second  order  error  bound 


(A'-A)^  >  0  , 
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CHAPTER  V 

A  direct  comparison  between  an  ordinary  first  order 
perturbation  expansion  for  a  quartlc  potential  and  a  perturba- 
tion In  the  path  Integral  formulation  for  the  same  potential 
will  be  made.   The  presence  of  certain  oscillating  terms  In 
the  path  Integral  formulation  will  be  eliminated  by  renormall- 
zlng  about  a  frequency  shifted  harmonic  term.   In  the  full 
N-body  problem  a  similar  technique  is  employed. 

A.   The  Ordinary  Anharmonic  Perturbation. 

The  one-body  velocity  auto-correlation  generator 

expressed  as  a  phase-space  integral  is 
1 


(5.1)   <e 


Rx(t)x(0) 


r      -  i-  Rx(t)x(0) 
>  =       e 


—  CO   _oo 


-BH(x(0),x(0)) 
e  dx(0)  dx(0)   ; 


In  this  equation  the  hamlltonian  is 


H(x(t),x(t))  =  Hq  +  H^ 


where 


1   '  2      1    2  2 
H^=  -^  mx  (t)  +  p-  ma)„x  (t)   and   H-, 


2(-2;^)x  (t) 


The  corresponding  equation  of  motion  is 

2 
(5.2)      mx(T)  +  mtOQx(T) ^  x^(t)  =  0 


with  solution  [25],  to  first  order  in  y,  given  by 
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x(t)    =    Xp^d)    +    yx-,(T),    where 


(5.3)    Xq(t)    =    (xQ(t)sin   WqT+   Xq(0)    sin   WQCt-T )  )/sin   oa^t 

and 

-,  -J  x^(t)-x^(0)cos    C0p,t    2 

(5.^)    x^(t)    =   ^   [x^(0)-    3Xq(0)(    -^ ^ ^)    ] 

sin   to„t 

_  X^Ct  )-Xq(0)coS     CJQt     o 

•[cos    Bu^qT-    cos    a)pT]+   ^   [( 3in   (.^t ^ 

2           XQ(t)-XQ(0)cos    tOgt        , 
-^3Xq(0)( 3^^   ^^t ^^4  ^^'^    ^''O'^-    ^^"^   '"O^^ 

subject  to  the  conditions 


x(0)  =  Xq(0)  ;     x^(0)  =  0 


x(0)  =  Xq(0)  ;     x^(0)  = 


=  0  . 


Therefore  (5.1)  becomes  after  expanding  to  first  order  in  u 


7   7   -  I  RxQ(t)xQ(0)   -|m6x^(0) 
e  e 


^0 


_00   —00 


-   imBa3QXQ(0)           ymBa)QxJ(O)         yRx^(t)xQ(0) 
■e  (1+   pTj ~   2 +...)D[ 


(5.5) 


^    (1+   2    ■•"    ^JM  ^0^^"^    3wQt+   j^  Wq    sin   w^t] 

24  0 


TcTe 


3  3 

Rwp.    sin   oj„t         9R   sin   to^t         3R   w^,    sin   w^t 


2ae 


2 


2aj„a    9 


Sa^e^ 


where    a   =   R   cos    Wp.t    +   m3         p 

r2    sin^^w^t 

^    =   '^0^"^^    -    4(R   cos    a3Qt    +   mB)^ 


B.   Path  Integral  Perturbation. 


We  now  ask  that  this  integral  give  the  same 
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result   as  the  path  expression 

-   i  Rx(t)x(0) 
<e  >  = 


-  I  Rx(t)x(0)  -  |m3x^(0) 
e  e 


—  oo      _oo 


(5.6) 


1  fir  2  2,^,   a|yx\o) 

2  mgCoaQX  (0)-  j2 ) 

e  exp[i 


5(t')(x(t')-x(0))dt' 


t  t' 

r 


+  i 


(coQ-f2^)C(t'  )x(t")  dt'  dt"  -  i 


t  t' 
r 


0  0 


5(t' )F(x(t"))dt'dt"] 


0  0 


t  t- 

1  rd?(t')  ]  rdx(t')  dx(0)  dx(t)  dx(0)  , 
0  0+ 


where  ^      represents  the  frequency  shift  term  about  which 
renormalizatlon  will  be  performed  and  F(x(t"))  contains  the 
cubic  force  and  linear  force  terms.   Simultaneously  putting 
the  kernel  into  the  quadratic  square  form  and  expanding 
one  obtains 


oo       oo 


(5.7) 


exp  - 


1 


-[ 


t  t' 
r 


_00       _00 


0 


2mBM  J   ?(t')?(t")  df  dt"-iRx(t) 

n   0.5     -  in3w^x^(0)   i  /  C(t ' )x(t' )dt ' 
-  1  R^x^(t)]e   2    0      ^   o 

P  t  t •  „  t  t' 

iaj„  /  /   C(t' )x(t")dt'dt"  -IQ'^    f   f    C(t' )x(t")dt'dt" 


C(f)dt' 


0 


o  0 


o  o 


2  H 
ymBcj^x  (0) 

(1+ i 

2H 


t  t 


-^  C(t' )x^(t")dt'dt"+ 


0  0 


tt 


to  within  a  Jacobian  transformation  factor  between  coordinate 
sets  (x(0),x(0))  and  (x(t),x(0))  which  is  essentially  the 
normalization  discussed  in  Chapter  II  and,  since  we  are 
going  to  renormalize,  is  not  included. 

See  Appendix  E  for  integration  of  linear  terms  in  the 
functional  integral  and  Chapter  VI  Section  C  for  further 
discussion  of  the  frequency  shift. 
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+  i 


t  t 


0  0 


f2  C(t  '  )x(t")  dt'  dt")  D[.  .  .] 


C.   Comparison  of  the  Integral  Forms 


We  will  now  compare  the  exponents  of  (5.1)  and 
(5.7)  without  the  frequency  shift,  although  the  shift  Is 
of  course   required  to  annul  the  secular  terms  In  general. 
Certainly  upon  x(0),  x(0)  Integration,  the  exponents  come 
down  to  the  Integrand,  and  hence  will  occur  together  with 
the  explicit  integrand  of  the  above.   The  required 
integrations  are  not  easy.   However,  our  major  interest 
is  in  long  time  dependence.   We  shall  show  that  in  this 
limit,  explicit  Integrands  become  Identical  without 
frequency  shift.   Hence  exponents  are  identical  too,  and 
further  integrations  are  not  necessary  to  establish 
identity  of  the  two  methods.   With  this  in  mind  expression 
(5.7)  can  be  written 

-  |(c^^x^(0)  +  2c^2^(°)^^^^'^^22^^^^^^ 
e 


00    oo 

f 


_00   _00 


(5.8) 


2  H 
ymBw„x  (0) 

(1  +  ^ 1 


t  t 


2 


r  f   w^y 


TW 


^  5  (f  )x-'(t")dt'dt" 
6   ^s      s 


0  0 


+  1 


t  t 

•  • 

J 

0  0 


fi^  K    (t')x  (t")  dt'  dt")  dx(0)  dx(t)  , 
s      s 


where  C  (f),  x  (t")   are  C(t'),  x(t")  evaluated  at  the 
s        s 

saddle  point  and  are  given  by  (2.9),  (2.10)  and  the c^j  are 
given  by  (2  .  l4  )  . 
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The  solution  of  equation  (5.2)  has  been  so 
constructed  that  its  zero-th  order  portion  is  exactly 
that  obtained  by  solving  the  variational  equations  of  the 
path  integral  for  the  one-body  oscillator  .   An  immediate 
observation  is  that  the  zero  subscript  can  be  placed  on 
all  X  or  X.   Thus  directly  comparing  the  kernels 


(5.9) 


12  2 

exp  -  2(c-j_-|_Xq(0)+  20^2^0  ^^^■^O^^''"''  ^22^0^^''^ 

12     2  2 

=  exp  -  p-  [{(BmcOj-,+  Rco„  cos  a)„t)/sin  bi^t} 

2       2        2        '2        2       2 
•(x„(0)+  x„(0)  cos   (ja„t  +  Xp^(O)  sin   w„t/a)„ 

1    2 

+  2Xq(0)Xq(0)  cos  W|-^t  sin  oj„t  /Wq)  -  2{(p-  Ra)„ 

12     2  2  2 

+  p-  Rw„  cos  oj„t  +  3m(jo„  cos  a)„t)/sin  w„t} 


(Xq(0)  cos  Wpt  +  Xq(0)Xq(0)  sin  oj^t  /^q  )  ] 


since 


(5.10)  Xq(t) 


T  =  0 


gives  x„(t)  =  x„(0)cos  w^t  +  x„(0)sin  to^t/w 


'0 


^0 


0 


^0 


^0^  -Q 


After  further  algebra,  this  can  be  reduced  to 


1    2  *  2 

exp  -  p-  RWq[Xq  (0  )x^  (0 )  cos  w^t  /  co^  sin  co  t 


+  Xq(0)  cos  cOgt/coQ  -  Xq  (0)Xq(0)/oJq  sin  oo^t] 

-  I  m3xQ(0)   -  I  m6coQXQ(0) 
•  e  e 


see  conditions  (5.^) 
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which  in  fact  corresponds  to 


(5.11)   exp  -  ^  R  XQ(t)xQ(0)  e 


-  ifn6xQ(0)   -  I  m3coQXQ(0) 
e 


since  Xq(t) 


(5.12) 


T-^^  gives 

XQ(t)     XQ(t)  cos  Oipt  -  Xq(0) 


00, 


sin  tOj-.t 


and  x„(t)  is  given  by  (5.10). 


It  remains  now  to  directly  compare  the  integrands; 


thus 


1  r,J 


3Xq(0)x^(0) 


-yRx^(t)xQ(0)  =  -yR{^[x^(0)xo(0) ^ ] 

o      ^Q  ^^ )  2  '2 

(5.13)-[-3Wq    sin    Bw^t+ooQ    sin   a)Qt]+    j—l — 5 3Xq  (0  )Xq  (0  )/aJn  ] 


^0- 


CJ 


0 


•l^r.    COS    3(jOQt    -    (jjp,    cos    tOgtlJ 


should  give  the  same  result  as 
t  t'  .2 

(5.1^)  -   ^    \    \ 
0  0 


-^  ?(t')x^(t")  dt'  dt" 


after  various  terms  are  renormalized .   Here  ?(t')  and  x(t") 
have  the  same  form  as  (2.9)  and  (2.10)  so  that  when  this 
integral  is  expanded  out  one  has  only  to  evaluate  the 
following  four  terms. 


t    t' 
f 


0    0 
1 


3co 


cos    uJq    (t-f)    sin^cogt"    dt '    dt" 


cos    oOgt    sin    4(jO|-|t         sin   tOgt      -, 
-[-    T-o": +  h         (■ 


0 


32a) 


0 


203 


0 


cos    2to„t 
2ooT         ■ 


sin  Wj-jt      -, 
8        ^^ 


cos    '^w^t  p      COS    (jo„t    sin    2(jj„t         sin   w^t 


^0 
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sin    co„t    cos    2aj„t  2    sin    a)„t 

0  ScOq 


t  t 


T 


cos    coQ(t-t')    sln^WQ(t-t")    dt'    dt" 


0    0 


1  1  13 

— PT  sin    4c0p,t-   — t"  sin    2wp|t    -   — ^rC  3    sin   aj„t    cos    w^t 


2    ^^''    --QO-         2    ^^''    "'"0"         ^    2^^    ""''    "^0"    """    "^0' 
9dcOq  doOq  9a)Q 

-   6    cos    (jo„t    sin   w„t] 

t    t  ' 

r     (  2 

cos    oaQ(t-t')    sin   tOgt"    sin   a)Q(t-t")    dt '    dt"    = 

0    0 

sin   a)„t                            ,            cos    2a)„t        sin   Wp^t    sin    2oo„t 
— T [COS    w^tiTi -   n )-   n 

3(jOq  0       ^(J^n.  40J„  4C0j-. 

I  COS    2a)j-,t  -,  cos    ^w„t 

-     cos     (J0r,t(-n -     n )-     COS     W^tC^'T -     ^o ) 

0       00)^  oWp  0       32co„  32Wj-, 

sin   w^t    sin    ^w^^t  cos    oa„t         cos    a)„t    sin    ^oj^t 

+   )1 !l-l    +   —[  -   - — 

32tOQ  3^Q  32tjQ 

,  cos    2a)„t  -,  cos    ^w„t  2    cos    tOj-jt 

+    sin   oj„t(-n —  -   n )    -    sin   oo^tCTT, —  -   ^c-^ )  ]+   5 

0    '8(j0q  8a3Q  0       32(j0q  32^0^  3^^ 

cos    to^t    sin    2co^t  .,  cos    2a)j^t  2    cos    w^^t    sin   w^t 

[ n +sin   aj„t(Ti -   — r. )]-   :t 

0  0  0  3(jJ„ 

t    t' 

2 
cos    (jOgCt-t')    sin   W|~,t"    sin   tOj-,(t-t")    dt '    dt" 

sin   oop^t      sin    ^oOp^t         cos    2oj„t  -,  cos    co^^t         sin   ^w„t 

=    _   !±_r    LL_    I  (         .       ^      _    .  -*-    '1 1      x_r_   !^ 

3oJq       ^         32a3Q         ^         5i;i^  ^^^  3a)Q      ^  32a)Q 

2  2 

sin   2a)p.t  sin   oj„t  sin   co„t    +2    cos    ca„t 

5 ^]  +  —   [+  ^^ 5 ^] 

2a)Q  ^0  3a)Q 

Therefore  equation  (5.1^)  after  algebraic  manipulation 
is 
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L 7- +      -,  p       ( COS    oj„t    -   x^(0)    sin   ojQtJJ 

.  T  .2  *  2  2 

Xq(0)         3Xq(0)Xq(0)    cos    w^t         3Xq(0)Xq(0)    cos    oogt 


OJ 


I 


ojp,    sm   co„t 


co„    sin   a)„t 


3  3 

x„(0)    cos    tOp.t  -,  cos    co„t    sin    ^co^t         sin   a3„t 

+    — o —){-  T^C \o ~  +         o,, 

•     3      ^-  3w„  32a)p,  003^ 

sm   Wj^t  -"    0  -^0  0 


sin   tOj-jt    cos    2(jj„t         sin   a)„t         sin   oo^t    cos    ^'^gt 


50) 


0 


32a) 


0 


32cj 


0 


P      cos    cOj^t    sin    2to„t         sin   ci)„t    cos    2co„t         sin   w^t         sin   w^t 


0 


0 


0 


0 


0 


2  ^ 

x^(0)x„(0)         x;?(0)    cos    o)„t  sin   a)„t    sin    ^o)„t         cos    2o)^t 

+    3(-^ o^ +   — o —)i o^^^[ — ^  + 


^0 


0)^    sin   o)^t 


3 
sin   o)„t 


30) 


0 


32o) 


0 


tO) 


0 


.     2 


cos    o)^t  sin    4o)^t         sin    2o)^t         sin   o)^t        sin   Qj„t 


IT 


0) 


0 


30) 


[- 


0 


0 


32o) 


0 


2o) 


-]- 


0 


0) 


-[- 


0 


3o) 


0 


2    cos^o)^t  Xn(0)x^(0)         2xn(0)x^(0)    cos    a)^t         x^(0)cos    a)^t 

o_])+  3(_o 0 ^  _^ g — _ o_  ^  _o _ u_) 


3co 


0 


2    •  ^ 

o)p,sin   o)„t 


o)j^sin   o)p,t 


3 
sin   o)Qt 


sin   o)„t    cos    o)„t         cos    o)   t    cos    2a)Qt         sin   bi^t    sin    2a)Qt         cos    o)Qt 


3o) 


0 


(0) 


0 


TT 


0) 


0 


"5 


0) 


0 


"5 


0) 


0 


cos  o)„t  COS  2o)Qt    cos  0)  t    cos  o)Qt  cos  '^w^t    sino)Qtsin4o)^t 


^0' 


iO) 


0 


32o) 


0 


32o) 


0 


32o) 


0 


^ 


cos  u)„t     COS  o)„t  sin  ^oj  t    sin  oj^t    sin  ui^t    cos  2o)Qt 


3u) 


0 


32o) 


0 


"B" 


0) 


0 


^ 


0) 


0 


sin  o)^t    sin  o)„t  cos  ^o)„t    cos  o)„t  sin  2o)„t    sin  oj„t 


0 


0 


0 


0 


0 


0 


32o) 


0 


32o) 


0 


2o) 


0 


2o) 


0 


;in  o)f.t  cos  2o)nt    2  sin   o)nt     ^n^^^        i 

0  ° ■ ^])+  ^-^ (-  ^^   sin  ^o)„t 


0 


'0 


sin  o)^t 


^0 


,3, 


sin  2o)p,t    sin  o)j-,t  cos  oi^t    2  cos  o)Qt  sin  o)Qt 
_ _  ;-  +  ^  ; 


6o) 


0 


3o) 


0 


30) 
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D.   Renormallzatlon  of  Terms 

We  now  search  expression  (5.1'+)  for  typical  terms 

which  correspond  to  those  of  (5.13).   The  first  such  term  is 

3xn(0)x^(0) 

(5.15)   -  {-^ )  (3  sin  3ajQt  -  sin  oj^t)  ; 

192c0q 

a  corresponding  term  arising  from  (5.1^)  is 
1    xp(0)x3(0)       ,   .   ,   , 

X  (0)x3(0)   Q  1  2  \ 

"  3(192)  J^ ^^  ^^"^   ^"^0^"  IT  ^^'^  "^0^"  ^^^   '^0^    ^^^   ^'^O^ 

since 

cos  (jO(-,t  sin  ojj-^t  =  2(l-sin  tOj,t)(sin  Sw^t-sin  cogt) 

3  1  3 

sin  co„t   =  -   jr   sin  S^Qt  +  jt-   sin  oa„t  . 

Therefore  the  terms  will  have  the  same  structure  provided  that 

3TTf2T  ^^0^°)^0^°^  sin^ojQt  sin  Si^^t 

can  be  balanced  by  the  secular  term  coming  from  the 
frequency  shifted  harmonic  contribution.   Another  contributor 
to  (5.15)  is 

1    ^q(°^^o^°^      .     •     u     . 

^  /  -I  n  o  \    cos  W^t  sm  HOJ„t 

2(192.;      co^  0        0 

which  can  be  converted  to 

(3)(4)(192)  ^  ^q(0)xI{0)    {sin  3a)Qt-  i  sin  a3Qt} 

leaving  also  a  portion  to  be  balanced  off.   Furthermore, 
other  such  terms  exist  and  can  be  handled  in  like  manner. 
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A  second  distinct  typical  term  is 
(5.16) 


3      Xq(0) 
~   Top    p~    (cos    3oj^t    -    cos    a3„t) 


a   correspondence      comes    from 


and   since 


11-4  cos    oj^t 

ETI92T  "2    ^0^°)    sin    a),t    ^^^"^   '^O^    ^°^    "^^^O^^ 
a)Q  0 


-   3 


cos   Wj-jt    cos    ^w„t    =   4('^os    S^i-st-cos   to„t ) 


+    (cos    cop^t    -    2    sin   co^t    cos    3co„t), 


we   obtain 


1      ^S(o) 


-    hdqo)    2~    ^*-^°^    3cOgt    -    cos    oigt) 

plus  a  portion  to  be  renormallzed .   The  remaining  two 
characteristic  terms  are 


(5.17) 

and 

(5.18) 


-1       Q       • 

JQ2-  Xq(0)Xq  (0)  (-3tOQ    sin    3coQt    +    Wq    sin   oa^t) 


Q         2  *  2 

j^  Xq(0)Xq(0)  (cos    3wQt    -    cos    OjQt)     , 


and  their  corresponding  components  are  treatable  in  analogous 

fashion.   Thus,  we  see  by  judiciously  balancing  terms  we  can 

get  the  integrands  to  be  identical. 

It  is  to  be  recalled  that  we  are  interested  in 

1 


<x(t)x(0)>  =  -  2  |p  <e   ^ 


Rx(t)x(0) 


R=0 


and  hence  can  discard  all  R  terms  in  the  examination  of  the 


-53- 


shifted  frequency  term  viz.  its  R  dependence  being  super- 
fluous following  the  differentiation  procedure.   Whereupon 
after  the  frequency  shift  term  of  (5.8)  is  integrated  and 
all  remaining  renormallzable  terms  are  collected  simple 

algebra  gives:        .^  .p 

2      2'       ^0^°^     3      3Xq(0)Xq(0)  2 

(5.19)  ^      =  ijoJqX„(0){( T—  sin  ajj-,t+  p cos  od    t    sin  a)p,t 

2       2 
3Xq(0)x„  (O)cos  cj^t  sinoagt 

+ +x^(0)cos-^u)^t)(-  ^-— - 

cos  a)„t  sin  ^^Qt    3  sin  w  t    sin  co^t  cos  2a)^t 
32wq  32a)Q  8to^ 

sin  w„t  cos  ^^p.^        cos  u^t  sin  2a)„t    sin  (jj  t  cos  w  t 

2 
3  sin  Wpt      Xq(0)x-,(0)  o  sin  Ugt 

-  5 ])  +3( sin  a)„t  +  x„(0)cos  aj„t)(-  — = 

^coq  coq        u     u       u      ^'^n 

sin  4w„t    cos  2(ji)-,t     -,     cos  a)„t     sin  ^oj^t    sin  2a3„t 

r U    ,   U        X  -|    U    r  U  ^1 

2  2  "2 

sin  a)„t    sin  oo^t    2  cos  w^t      x„(0)x„(0)     ^ 

-   7-, L-  — 5 Tr^ ^^^  3( 0 sm  a)„t 

2 
2Xq  (0  )Xp,  (0  )cos  w^t  sin  w^t     -,       P      sin  oj„t  cos  co„t 

+  -2 +  Xq(0)cos  ^ot)(   3     [  i^ 

0  u       u 

cos  Wpt  cos  2a)Qt    sin  w^t  sin  2a)Qt      cos  w  t 


^COp,  IjTI)" 


0  ^'"O  "^0 

cos    a)p.t    cos    2to„t         cos    w^t        cos    oo^t    cos    ^w„t        sin   aj„t    sin   4a)j-,t 
■^  ffZ;;^  32^  "^  32^;^^  "^  32^;;j^  ^ 

cos    fcOj^t  COS    a)„t    sin   ^w„t         sin   oo^t         sin   oo„t    cos    2co„t 

3co^  32c0q  Bw^  BIo^ 
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sin  cOpt         sin   w   t    cos    '^'^^t         cos    cj„t    sin    2oj„t         sin   Wj^t 
32^  +  321;^^  +  2^^  +         20)0 

sin   co„t    cos    2co„t         2    sin   oa„t  .^  sin    4aj„t 

-^, — s^"^  '<^°)(-  ^^ 

3 
sin    2coQt         sin   co„t    cos    co^t         2    cos    oj„t    sin   w„t 

-  2        ~  2  "^  2 

6ajQ  3a)Q  3wq 

P                        sin  oOp^t    t  cos  00p.t 
/{(xQ(t)-XQ(0)xQ(t)  cos  cjot)(   2^  " 2 ^) 

P             a.    t  COS  oo„t    sin  2oo„t   sin  ajp.t 
+  (Xo(0)Xo(t)-x2(0)  cos  C0ot)(| ^^  +  -.^^^^-^ ^)} 


Now  as  t  ->■  0°  the  numerator  has  oscillating  sine  and  cosine 

terms  whereas  in  the  denominator  the  secular  t  term  goes 

2  t 

to  infinity  indicating  that  Q.      ->  0  in  the  long  time  limit. 

We  explicitly  write 

2  1 

(XpCt  )-XQ(0)xQ(t  )cos  ajQt)(-  J   t  cos  w^t )  +  (x^COxQCt) 

2  1 

-  x„(0)  cos  oo„t)(p-  t)(l  -  cos  to„t) 

2 
=  [x^  (0  )Xo  (0  )  (sin  w„t-  sin  u)„t  cos  co„t  )-x„  (0  )x„  (0  )sin  w„t  cos  aj„t 

'22  t 

-  x„(.0)  sin   (jj-,t  cos  co„t  /a)„] 


since  x„(t)  =  x„(0)  cos  oa^t  +  x„(0)  sin  w„t  /co^  • 
The  coefficient  of  the  t  term  is 

2 

x^CO)  sin  oj  tCx^  (0 )  (1-cos  <^Qt)-  Xq(0)  cos   (jOgt 


(5.20) 


-  x„(0)  sin  (jo„t  cos  (ja„t  /Wj-,]  , 


t 
this  can  be  seen  to  be  true  for  the  imaginary  time  case 

as  well  viz.  see  Appendix  D.III. 
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and  the  sin  a)„t  factor  can  be  cancelled  since  each  term  in 
the  numerator  has  a  like  factor.   Therefore,  we  must  now 
check  the  coefficient  to  make  sure  that  it  does  not  go 
through  zero. 

E.   Examination  of  Secular  Coefficient. 

If  the  coefficient  did  become  zero  for  some  values 
of  x„(0)  and  Xj,(0)  we  would  have 

2  ^0^°^  ? 2 

(5  .  21  )x„  (0 )  (1-cos   oJnt-    cos    a)p,t)-   cos    w^t    /l-cos    (jo„t    =    0: 

U  U  U  (jO„  u  u 

transposing,  squaring  and  collecting  like  terms  gives 

(5.22)   c'^  +  2q  c^  -  c^  -  2q  c  +  q   =   0 

2       2      *  2      2 
where   c  =  cos  co-^t   and   q  =  x„  (0 )/ (x„  (0 )+  x  (0)/a)„)  . 

Upon  examining  equation  (5.22)  one  finds  that 

the  roots  of  the  quartlc  are  complex  and  hence  the 

coefficient  of  the  t  term  cannot  reach  zero  for  any  real 

values  of  cos  w^t.   Thus  the  two  perturbation  forms  give 

the  same  result. 


t 
See  Appendix  F, 


-56- 


CHAPTER  VI 

The  path  Integral  representation  Is  now  called  upon 
to  handle  the  N-body  problem  of  transport  in  a  Debye-type 
crystal  with  nearest  neighbor  cubic  forces.   To  guarantee 
normalization,  start  with  a  system  for  which  an  exact 
solution  is  known  and  turn  on  the  remaining  potential  from 
-°°  <  t  <  0  as  lim  $(t)  =  $„  +  $„  e^   (where  e.g.  $„  = 
harmonic  and  $   =  quartic).   The  density  subject  to 

y 

3p/3t  +  v9p/9x  +  v8p/3v  =  0  will  have  its  normalization 
maintained  since  8/3t    p  dx  dv  =  0  and  the  initial  entropy 
is  unchanged.   Then  at  time  =  t  one  inquires  as  to  the 
correlation  of  the  one-body  phase-space  density  with  that 
at  time  =  0,  allowing  the  calculation  of  transport  coefficients 

In  particular,  we  will  determine  an  inequality  for  the 
velocity-velocity  correlation  function  in  the  infinite 
imaginary  time  limit.  For  nearest  neighbor  type  forces 

(6.1)  f  (x(t"))  =  -  y  A  .(x  -x.)^  +  y  m  fi^.(x  -X.) 

s''''      ^sj'sj'     ^     SJSJ 

where  x.  -   x.(t");  s  =  j  +  1  and  the  corresponding 

Hamiltonian  is 

p2 

(6.2)  E-l^.      ^_1a   (x^-x)^-   Ijm^^x^     2^ 

s        s,j       ^       ^       S,J 
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A.   Normal  Mode  Transformations 


t 


To  transform  to  normal  mode  variables  one  writes  the 
deviation  from  particle  position  as 

(6.3)     Xj.(t)  =   I   -^exp   ^^^   x(k,t)  +  j 

and  time  differentiating  gives 


(6.4) 


;.(t)  =  ):   1  exp  ^%ii^  x(k,t)  . 


M 


J       k 
Likewise  the  conjugate  variable  is  expressed  as 

(6.5)     C,(t)  =  I  7^  exp  ^^      (k,t) 


J      k 
with  lattice  modes  j  =  1,2,...,M  and  k  =  1,2,,..,M, 


To  apply  these  transformations  use  will  be  made  of 

,   2Trljk/M   ^   -27Tij£/M 
(6.6)     y  ^  e        -7^  e         =6, 


and 
(6.7) 


4  7m 
J 


7n 


ki 


A,  „      2iTln£/iyi      -2TTlnk/M  ^ 

V        kX/  f^ 

I      —nr-  e  e  =   mco 


k,Jl 


M 


n 


We  will  again  want  the  Integral  from  t  =  0  to  t  =  1' 


lim  [ 
a)->0+ 


X   e      .  .  .  > 


X(-B+l-ln  p)p 


X(B-B+l)e  ^. . .] 


with   X  =  e^-^  exp  -  ^  ^  r[x.(t)  -  x  (0)]^. 

Thus,  applying  the  transformation  to  the  Gaussian 
velocity  distribution 


t 


See  Appendix  A 
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(6.8) 


-6(j)(x(0)) 


=    e 


-    I   ^  m  xUo) 


(6.9)    =   exp   -   I   I  m3[I    7^  ^^^±3^/^   ;,^^^o)f 


J 

1 


k 


exp    -    p-  m3      I 


j,k,k' 


1      2.1Jk/M    ^2TriJkVM    •(i,^o);(k',0) 


=    exp    -   ^  m3       I         6,       ,   ,    x  (k  ,0  )x  (k  '  ,0  ) 
k,k'      ^'"^ 


(6.10)    =   exp   -   4  ^   m6x(k,0)x(-k,0) 

k 

B   for   the    path    is    6<})(x(0))    +    B(J)(x(0)) 
t 


-   1 


1 
m 


I    5.(t')(x    (f)    -   X    (0))    dt' 

J  J  J 


D    J 

0    t' 
f 


2^ 

m 


t       0 
f 


0    -°° 


^     5     (f)     f     (x(t"))     dt'     dt" 
S  ^  ^ 


I    C,(t')     f     (x(t"))     dt'     dt" 
s        ^  ^ 


(6.11)      -  i 
m 


t    t 


0    0 


I    i    (f)    f    (x(t"))    dt'    dt" 
s      ^  ^ 


for   t '    <    0 


for   t'    >    0 


for   t '    >    0 


where    Its    potential    Is    giver,  by 


(6.12)      <f(x)    =      I      1   I^.(x^-x^.)^   -      I      i  mJi2^.(x^-xJ^ 


SJ         S        J 


SJ         S        J 


The  purely  quartic  part  in  the  mode  variables  will  be  turned 

Yt 

on  slowly  from  t  =  -°°  to  0  via  e'   while  the  harmonic  part 
of  the  quartic  Interaction  arising  from  the  normal  mode 
transformation  will  be  lumped  with  the  ordinary  harmonic  term, 

Hence,  applying  (6.3)-(6.5)  one  may  write  B  in  obvious 
notation  as 
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(6.13)    Um^ll   4e2-iJ^/^  x(k,0)]2 


k 


^    6eVt      J      1    ~^^^[    ^    (e2-i=K/"-e2''lJ>^/")Si^.(s-J)]'' 


S5J 


1    -       V         n2    rV^    2TTlsk/M      2TTljk/Mvx(k,0)    ^,       ...2 
-   2    ^      ^      mf^      [^(e  '     -e        ^    '     )      y^        +(s-j)] 


-^       4\i  7m 
0  J  k 


-x(k,0))    dt' 


+   i 

m 


0      t 


.00       .00 


Id 

s   k 


1      27Tlsk/K 
7m  ^ 


C(k,t'))(I    A_.[^    (e 


27Tlsk'/M 


^       sj 


J       --    k- 

_     ^2TTljkVM) 


.    £(^l.i::l.(3-j)]3_    ^   ^nln    (e2-skVM_^2.1JkVM),,(kl^ 


J         ^J    k' 


7m 


+    (s-j)])    dt '    dt"    +   ^ 


m 


0    -0° 


0       0 


while   B   is    given  by 

rt^    n;n    1    V      orV       1      2Trijk/M    '    ,     ^x-i2_^    1  r       «       o^Vr    2irisk/M 

(6.14)    ^  I   m^ll    j^  e        ^           x(k,0)]    +    ^  A      A      6(Ue 

j          k  s,j        ^      k          ^_ 

t 

2Trijk/M-,x(k,0)    ^.       .vv2       .     f    v  ^  V  1      2TTijk/M    ^,,     ,  ,  >,  ^ 

-   e                 ^-Tm —  +(s-j))    -   1        Hz   7m   e  ?(k,t')) 


J  k 


t  t 


k '  J  •'   s  .  .1  k 


0  0 


.^2^isk/M  ^(k,t'))(A^.  I      -7I  e^^iJl^'/M  ^d^.^t"))  dt  '  dt" 
The  simplification  of  B-B+1  proceeds  as  follows.   The  term 


(6.15)   I  B  I      mQ^ 


s,j    ^J  k 


(^^^2uisk/M_  ^2.iJk/M3  x(k^  ^  (3_.))2 
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reduces    to 
(6.16) 

(where 


I    mB    e 
k 


+  2TTik/M 


!;2^x(k,0)x(-k,0) 


;2 


s,J 


mfi^^.      2iTlkj/M      -27Tiks/M 


I      —^   e        e         =  mQ"-)  , 


and  a  linear  term  of  the  form 


(6.17) 


2TTiks/M 


i  I   m6/M   Q^    e  x(k,0) 


k 


together   with    Its    symmetric    counterpart    in   -k 


Skipping    for   a  moment    the    thii'd   and    fourth   power   terms, 
t    t' 


(6.18)    - 


m 


0    0 


1,(1   7^  e 
s,j    k 


2TTlsk/M 


C(k,t'))(A3.    l^    ^e 


2TTijk'/M 


x(k' ,t"))    dt'    dt" 


becomes 


(6.19)   -i 


t      t' 

r 


A    .      2TTisk/M      2TTijkVM 


0    0      f 'J 
k  ,k ' 

t    t' 


M 


5(k,t' )x(k' ,t")dt'dt" 


(6.20)    =   -±    I    (    \   a)^^(k,t')x(-k,t")    dt'    dt" 
k   J    I 


I.e. 


A 


^  .      M 
The  third  order  term 


0    0 


sj       27Tikj/M      -2Triks/M  2 


=  mco 


In  demonstrating  that  the  frequency  shift  goes  to  zero 
finite  terms  such  as  these  are  found  in  the  denominator 
of  the  expression  for  9.      (see  Appendix  D.III  or  equation 
(5.19))  and  as  time  becomes  infinite  are,  in  comparison 
with  the  secular  terms,  of  no  consequence. 
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A    .  2TTisk/M  2iTisk'/M      27TlJkVM 

(6.21)         I        -^  C(k,t')    e  il    [e  -e  ] 

s,j,k  k- 

__        ^        ^, (,,,.,    e''''''^{I    [e2-i^^VM_^2.iJkVM3x(k^^ 
s,j,k  k' 

+    (s-J)}{I    te2--l^"/M_^2.1JkVM3    x(k%t::i  ^    (^_.)j 


k" 


^^    ^^2.iskMVM_^2.1Jk'VM^    x(k"'    t")    ^    (^_.)j) 


k'" 
simplifies    to 


(6.22)  ^  e^^"   -^  C(k,t')x(k' ,t")x(k",t")x(k"',t")e^^^^^/'^ 
s,j,k  M 

k'  ,k",k"' 

2u^(sk'+sk"+sk'"  )  2TTi(jk'+sk"+sk"'  ) 

•  {e  -   e 

27ri(sk'+jk"+sk't' )  2ui(  jk '+jk"+sk  "'  ) 

-e  +   e 

27T^(sk'+sk"+jk'i'  )  27Ti(jk'+sk"+jk"'  ) 

-e  +   e 

2TT^(sk'+Jk"  +  jk"')  2T7|(jk'+jk"  +  jk'") 

+  e  -   e  } 

Once  again  odd  powers  In  x  are  discarded  (recall  ?  contains 
X  explicitly)  and  second  power  even  terms  will  be  lumped 
with  the  harmonic  contribution.   Hence,  In  like  respect  the 

fourth  order  term 

,    .      2ulsk/M    27TlJk/M    n  n^         ;i 

(6.23)  I     I   B  A  .[^(e       -  e       )  ^^^^  +  (s-j)]^ 
s,j       ^  k 

Is  expressible  as 
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A    . 


(6.24)        I  e^^   J  3   -^  x(k',0)x(k",0)x(k'"  ,0)x(k""  ,0) 

k",k'"  ,k"" 


{e 


27Ti-(sk'  +  sk"  +  sk'"+sk"" 


M 


2TT^(sk'+jk"  +  sk'"+sk''" 


-e 


M 


2TTi-(sk'+sk"+jk'>'+sk"" 


-e 


M 


2TT^(sk'+jk"+jk'"+sk"" 


+e 


'M 


2^^(sk'+sk"  +  sk"'+jk"" 


-e 


M 


2TT^(sk'+jk"  +  sk'"+jk"" 


+e 


'M 


2TT^(jk'+sk"+sk'"+sk''" 


M 


+   e 


+    e 


2u^(jk'+jk"+sk'i'+sk'"i 
2^i(jk'+sk"+jk"'  +  sk"" 


-   e 


2wi(jk'+jk"+jk'"+sk"" 


+    e 


2^^(jk'  +  sk"  +  sk"'+jk"" 


2TT^(jk'+jk"  +  sk'"+jk""   ) 


M 


2TTksk'  +  sk"+jk"'+jk'<"   )        2ui(jk'+sk"+jk'"+jk""    ) 


4e 


M 


-   e 


M 


2ui(sk'+jk"+jk'"+Jk""  )        2^i(jk'+jk"+jk'"+Jk""  ) 
-e  +   e  } 


The   kernel 


t   t 


(6.25)    exp    -    2^[ 


■.    L    i    k  k ' 


ooJ  ^ 


1      2TTlJk/M 


?(k',t"))dt'dt"-ir    J     m    ^  e^-J-/"'    x(k,t))(I^    ^ 

^2.iJkVM    ^^,,^,.)    ,,,    _    1    ^2    ^^^    _!    ^2.1Jk/M   ;^,^t))(I      ^ 

j    k  k' 


.^27TlJkVM   ;(i,,^t)]    exp    -   I      I      B   A    .(ICe 

s,j  '^    k 

t 


2TTlsk/M_g2TTijk/M-| 


M 


^^^^^)'expl    f    ^(^^e2^^J^/"c(k,t.)) 


J    k 
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t  t' 
f   f 


(I      T^e^^-J'''/"      i(k-,t.)dt.)    expM  in   7! 


7m 

2TTlsk/M 


0    0 


s,j    k 


C(k,t'))(A^.    I^    ^e 


1      2TTijk'/M 


x(k',t"))dt'    dt" 


becomes 


t    t 


(^•26)    -P    -    2^   t 


0    0 


I        6^   _^,    5(k,t' )5(k',t")dt'    dt" 


k,k 


t 
f 


-    ir 


I        6,       ,   ,    x(k,t)C(k',t')dt' 
k,k'      ^' 

^        6      _        x(k,t)x(k',t')] 


0 

1      2 

^"^      k^k-    "^'-^ 
+  2TTlk/M      2 


-m^l   e-''"^'^/'''   oo^x(k,0)x(-k,0)      i    /      ^      ^  (k  ,t  •  )x  (k  '  ,t '  )dt 


•e 


k 


0    k,k' 


t   t 


2 


•e 


i/m   f   f      I   mu      C(k,t' )x(-k,t")dt'dt" 
0    0k 


or   finally 


(6.27)    exp   -  ^  n 


2m6 


t    t 


0    0 


-   ir 


C(k,t')?(-k,t")    dt'    dt" 


1      2' 


k 


x(k,t)C(-k,t')    dt'    -  J  r   x(k,t)x(-k,t)] 
aj^x(k,0)x(-k,0) 


qV      +2uik/M      2 


i    );    /   C(k,t')x(-k,t'  )    dt' 
k   0 

t   t'      p 
±    I    f   f     u)^(k,t' )x(-k,t")dt'dt" 
k   0   0 


where    the   kernel    is    equally    symmetric    in   k   and   -k , 
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B,      Establishment    of   Inequalities 


The    full    inequality    is 


(6.28) 


A    e-B   > 


X(B-B+1)    e 


-B 


(6.29)    = 


oo  oo 

f 


—  00  _00 


exp  -  75-^   y  [ 

k 


t    t 


0    0 


?(k,t' )5(-k,t")    df    dt" 


t 
r 


-   ir 


1      2- 


x(k,t)?(-k,t')dt'    -   J  r   x(k,t)x(-k,t)] 
-mg)]    e- 


+  2TTik/M      2    ,.     ns     /    ,     ^N 
u)   x(k,0)x(-k,0) 


t    t 


i    ^    /    C(l^,t' )x(-k,t  •  )dt'    ^l    f    f      mto    5(k,t' )x(-k,t")dt'dt" 


k    0 


k    0    0 


(1+   I      1      A    -BEI    (e 
s,j      ^J      k 


2TTisk/M        2TTi,ik/M^x(k,0) 


-   e 


Vm 


+    (s-j)]' 


.Nn2 


s,j  '^    k 


1 


.  Nn2 


±      V        nn2    rVr    2TTisk/M        2TTijk/M,x(k,0)      ^    ,       ... 
-    2      A      m3f2      [^(e  -^    -   e        ^         ) — ^.j^     +    (s-j)] 

s,j  ^    k 


1 

m 


t    t 


0    0 


r    /v      1      27Tisk/M 
s,j    k 


C(k,t'))(A^.    l^    ^e 


1      2    ijk'/M 


x(k' ,t"))dt 'dt"-   - 


0      t' 

yt" 


m 


_oo     _oo 


V^V      1      2TTisk/M 
s    k 


C(k,t')) 


(I    A    .[    I    (e2-3kVM_    ,2.iJkVM)XOcVLl).^(3_.^33    _ 


^      sj 


k 
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-   I   mf^cn-LZ  (e      '  -  e    "^      '     )  /' '-   +  (s-j)])dt'dt' 

.1    ^^    k'  "^^^ 


1 
m 


t   0 


0  -0° 


m 


t  t  '  . 

r  f  .  ^ . 

.  .  .  TT  d?(t  '  )  dx(t'  )  . 

0  0  ° 


t 

After  a  good  deal  of  tedious  algebra   the  inequality  can  be 
written 


(6.3C)       I    Xe   ^. . .    > 


X(B-B+l)e 


-B 


7 


—  00  _co 


00 

f   exp    I-   \  I    (c^^(k)x^(k,0)+   2c^2 ^^)^^^'° ^^ ^^'^ ) 


2 


+  2iTik/M      p      p 
+   C22(k)x''(k,t))](l+    I   mBe  (w  +f2    )x(k,0)x(-k,0) 

k 


I  i  3'|'(s,j,k,k',k",k'")    x(k,0)x(k',0)x(k",0)x(k'»',0) 


s,j  ,k 
k' ,k",k' ' 


+  ^    [m6(x(k,t)x(-k,t)-x(k,0)x(-k,t)cos    cot  );|j(s  ,j  ,k,k  •  ,k"  ,k '"  ) 

k 


(b^(t)  +  l)+  m3(x(-k,0)x(k,t)-x(-k,0)x(k,0)cos   tot) 


'ij;(s,j  ,k,k' ,k",k'"  )(bp(t)  +  l)+   ^  r  (x(-k  ,t  )x(k,t  )cos    cot 


x(k,0)x(-k,t))(b    (t))+   4  r(x(-k,0)x(k,t)cos    wt 


-   x(-k,0)x(k,0))(b.(t))]    +  I  e  (s,j,k,k',k",k'") 

s,J,k,k',k",k'" 


{[x(k,t)x(-k',t)x(-k",t)x(-k''',t)][a^(t)]   + 


t 


See   Appendix   D 


-66- 


+    [x(k,t)x(-k'  ,0)x(-k",0)x(-k'",0)][a5(t)]  +  [x(k,t)x(-k',0) 


■x(-k",t)x(-k"'  ,t)+   x(l^t)x(-k'  ,t)x(-k",0)x(-k'"  ,t) 


+x(k,t)x(-k' ,t)x(-k",t)x(-k'" ,0)^la    {t)^+    [x(k,t )x(-k ' ,0 ) 


x(-k",0)x(-k"' ,t)+   x(k,t)x(-k' ,0)x(-k",t)x(-k'" ,0) 


+  x(k,t)x(-k' ,t)x(-k",0)x(-k"' ,0)][a^(t)]-    [x (k  ,0 )x(-k ' ,t ) 


x(-k",t)x(-k'",t)][a^(t)]-    [x(k,0)x(-k',0)x(-k",0)x(-k'",0)] 


[ag(t)]-    [x(k,0)x(-k' ,0)x(-k",t)x(-k'" ,t)+x(k,0)x(-k' ,t) 


'x(-k",0)x(-k'" ,t)+    x(k,0)x(-k' ,t )x (-k" , t )x (-k '" ,0)][a^(t)] 


-[x(k,0)x(-k' ,0)x(-k",0)x(-k'" ,t)+   x(k,0)x(-k' ,0)x(-k",t) 


•x(-k'" ,0)+   x(k,0)x(-k' ,t)x(-k",0)x(-k'" ,0)][aQ(t)]} 

P      +27Tik/M 
-    I   Q.    {e  -l)[{m3(x(k,t)x(-k,t)-   x  (k,0  )x  (-k  ,t  )cos    tot 


+   J  r(x(k,t)x(-k,t)cos   wt   -   x(k,0  )x(-k,t ) }  {c-j_(t )  } 


+    {m3(x(k,t)x(-k,0)-   x (k,0 )x(-k,0 )cos    wt)+   \  r(x(k,t) 


x(-k,0)    cos    ojt    -   x(k,0)x(-k,0)}{cp(t)}]    D[  .  .  .  ]    ; 
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where  in  the  kernel  the  summation  over  k  stands  for  k,k',k", 
k '"  terms  and  plus  or  minus  values;  moreover,  D[ . . . ]  means 
in  each  summation  group   one  integrates  over  the  appropriate 
dx(k,0)  or  dx(k",t)  whatever  the  case  may  be.   Such  terms 
as  c^^(k),  c^^^^)'  C22(k),  4i(s  ,  j  ,k,k  '  ,k"  ,k '"  ) ,  b^(t),  b2(t), 
b^Ct),  bi^(t),   e(s,j,k,k',k",k"'),  a^(t),  a2(t),  a^Ct),  a|^(t), 
ai-(t),  ag(t),  arj(t)  ,    ag(t),  c^(t)  and  C2(t)  will  be 
identified,  but  first  some  of  the  terms  leading  to  this 
inequality  are  clarified.   The  combined  three  third  order 
terms  including  the  phase 


(6.31) 


m 


dt 


dt 


"   e''*"    III   ^i  e^-^=^/" 

s   k 


C(k,t')) 


a   IaI    (e2-skVM_  ,2uiJkVM)2LQ<l^,(3_j)]3 
j        J    k'  "^ 


m 


0 
f 


2^ 

m 


0    -^ 
are   reducible   to 


t   t 

•  r 

0    0 


(6.32) 


s,j,k        M 
k'  ,k",k"' 


A    .        27Ti^(k+k'+k"+k"' )      2^^(  jk '+sk"+sk  "•+sk) 
-^  {e  - 


-e 


27T^(sk+sk'+jk"+sk'*'  )        2u^(sk+jk '+jk"  +  sk '"  ) 
-e  +   e 


M 


2TTi(sk+sk'  +  sk"+jk'i'  )        2^^(sk+jk '+sk"+jk '"  ) 
-e  +    e 

2Tri(sk+sk'+jk"+jk'"  )        2TTi(sk+Jk '+jk"+jk '"  ) 
+e  -   e  } 
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rg  (x(k,t  )-x(k,0)cos    cot)  r    ^x(k,t)cos    ujt-x  (k  ,0  )  ^ -^ 

sin   cot  sin   cot 

•{[x(-k' ,t)x(-k",t)x(-k"' ,t)][-   I  -  cos    cot+   ||-  sin   cot 

_   I  sin   cot    +   ^  sin   a)t    +    sin   4(.+  cos    a)t    ^    5    sinco^t    cos    2a^t 

_    sin   oit^cos    ^cot    _    cos    a)t    sin    2a)t-j^    [x  (-k  '  ,0  )x  (-k"  ,0  )x  (-k 'm  ,0  )  ] 

r3  fJJ    cos    ojt    ,    sin    ^cot         sin    2ajt         1.3,                ^.,    sin   2oot 
•L-g-  :^ +   gg g 3    sm   a)t    cos    ajt+   ^ 

sin    4cot         sin    2cot         sin   cot    cos    oot         sin   2cot-, 
-^  ^  _  +  _         J 

+  [x(-k'  ,0)x(-k",t)x(-k"'  ,t)+x(-k'  ,t)x(-k",0)x(-k"'  ,t) 

2 
+   x(-k',t)x(-k",t)x(-k'",0)]    [_   11    sin   2    t    +   ^  "^   ^^"   ''^ 

2                                                                                                 2 
,     3    w   cos    cot    ,    1       .       n    .     ,    sin    2cot    cos    ^cot         sin   ojt    sin   2cot 
+   -s-  +   r=r  sm    2cot    +  


F  y "    3    ^^'^    ^'^^    "   9^ 1 


o 


2  2  2 

sin   cot    sin    4cot         cos    cot    sin    ^cot         cos    oot    sin    2cot         sin    ^cot 
+  g-g  g^  +  g  Y^ 

-   ^  sin    2cot]+    [x(-k' ,0)x(-k",0)x(-k"' ,t)  +  x(-k' ,0)x(-k",t) 

.    x(-k"',0)+x(-k',t)x(-k",0)x(-k'>',0)][^l'^   ^^^g°^    ^^^ 

2 
sin   cot    cos    ^oot         3    co   cos    cot         cos    cot    sin    4cot         cos    cot    sin    2cot 

9^  F  y  9^  5 

2 
sin   oot         sin   cot    cos    (ot         sin   cot         sin   cot    cos    2cot 


3         -  3  12  12 

sin   cot    sin    4cot    ,    cos    cot    sin    4cot    ,    cos    cot    sin    2cot 
9^ +  95  +  5 

3  2 

sin   (ot         2    sin   cot    cos    oot -i  ■, 

+  ^        +  3  ^-^ 

which  with   reciefining   the   phase    anci   simplifying  becomes 
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(6.33)  I  9(s,j,k,k',k",k"'){[x(k,t)x(-k',t)x(-k",t) 

s,J  ,k 
k'  ,k",k"' 

•x(-k'",t)][a^(t)]+  [x(k,t)x(-k',0)x(-k",0)x(-k"',0)] 
•[a2(t)]+  [x(k,t)x(-k',0)x(-k",t)x(-k"',t)+x(k,t)x(-k',t) 
•x(-k",0)x(-k"'  ,t)+  x(k,t)x(-k'  ,t  )x(-k"  ,t  )x  (-k '"  ,0)][a2(t)] 
+[x(k,t)x(-k' ,0)x(-k",0)x(-k'" ,t)+  x(k,t)x(-k' ,0)x(-k",t) 
.x(-k'",0)+  x(k,t)x(-k',t)x(-k",0)x(-k'",0)][a^(t)] 
-[x(k,0)x(-k' ,t)x(-k",t)x(-k'",t)][a^(t)]-  [x(k,0)x(-k'  ,0) 
•x(-k",0)x(-k"',0)][ag(t)]-  [x(k,0)x(-k',0)x(-k",t) 
•x(-k'",t)+  x(k,0)x(-k' ,t)x(-k",0)x(-k'"  ,t)+  x(k,0) 
•x(-k' ,t)x(-k",t)x(-k"',0)][a^(t)]-  [x(k,0)x(-k',0) 
•x(-k",0)x(-k'"  ,t)+  x(k,0)x(-k'  ,0)x(-k",t)x(-k'*'  ,0) 
+  x(k,0)x(-k',t)x(-k",0)x(-k'" ,0)][ag(t)]}. 

Likewise  the  fourth  order  term  may  be  rewritten  as 
well,  and  since  all  summation  indices  are  merely  dummy, 
one  may  let  k'  ^  k;   k"  ^  k';   k"'  ^  k" ;   k""  -^   k  •" 
which  makes  this  term  notationally  compatible  with  the  third 
order  terms. 

.   1    ~     .    2uisk/M    2uijk/M   „.,^  o^         u 

(6.34)  I      ^  3  A   [  ^  (e       -e       )  ^Li^  +  (s-j  )  ]^ 
s,j       'J   k 

(6.35)  =   I    ^  3i|j(s,j,k,k',k",k'")x(k,0)x(k',0)x(k",0) 

5,3  ,k 
k',k",k'"  .x(k"',0) 

where 
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..  s 


(6.35)'      e(s,j,k,k',k",k 


A    .         2TTi^(k+k'+k"+k"' ) 


'")    -    -^   {e 


"M 


2Tr^(jk'  +  sk"  +  sk"'+sk 


-e 


2TTi(sk+jk'+jk"  +  sk"' 


+e 


M 


2TT^(sk+jk'+sk"  +  jk'" 


+e 


M 


-e 


2^i(sk+jk'+jk"+jk'" 


-e 


ij;(s,j,k,k',k",k 
2ui-(jk+sk'+sk"  +  sk'" 


+e 


+e 


2ui-(jk+jk'  +  sk"+sk'" 


2iTi-(jk+sk'+jk"+sk'" 


-e 


+e 


-e 


-e 


2TT^(jk  +  jk'+jk"  +  sk'" 

2TTg(jk+sk'  +  sk"  +  jk'" 
2^^(jk+jk'+sk"+jk'" 

27Ti(jk+Sk'+jk"  +  jk'" 


+e 


27Ti(jk  +  jk'+jk"  +  jk'l' 


2TTi.(sk+sk'+jk"  +  sk'"  ) 


-   e 


M 


-    e 


2TTi-(sk+sk'+sk"+jk'"  ) 


+   e 


2TT^(sk+sk'+jk"+jk'") 


A    .        2TTi|(k+k'+k"+k'"  ) 


")    =   ^   { 


-M 


M 


2 


2TTi(sk+jk'+sk"  +  sk'" 


-   e 


M 


-    e 


2ui-(sk+sk'+jk"  +  sk'" 


+   e 


-   e 


+   e 


2u^(sk+jk'+jk"+sk'" 
2TTi-(sk+sk'+sk"+jk'" 
2^i(sk+Jk'  +  sk"+jk'<' 


+    e 


2w^(sk+sk'+jk"+jk'" 


-   e 


2^i(sk+jk'+jk"+jk'" 


}; 


a^(t) 


e.^{t) 


g         .3    to   cos    bit  sin    ^(iot    cos    a)t 
sm   (jot                   ' 

5    sin   oat    cos    2a)t  sin   cot    cos    ^ojt         cos    cot    sin    2a)t  v 

+   2l|  9^                                     5                   ' 

2 

°         ,3co   cos    cot  sin    2cot  >, 

C ^             +  3—^ 


sm   cot 
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a^Ct) 


sin    2ojt         sin    4a)t         3    (jJ   cos    tot 


r=0  .    5    ^ 

sm   wt 


TF 


T^  F 


Y 


sin   2tot    cos    4(jot         cos    2ajt    sin    4tot         sin   cot    sin    2a)t 


96 


"9F 


12 


,    cos    cot    sin    2cot  v 
+   ^ ) 


a^(t) 


r=0  .    5   4. 

sm   cot 


B         /      sin   cot         sin   cot    cos    2cot         sin   cot    cos    ^(ot 


12 


~2T 


"9F" 


sin   cot    sin    ^cot         3co    cos    cot    ,    cos    cot    sin    2cot    ,    cos    cot    sin    4cot>, 
9^ 8 ^  3 9^ ^ 


a^Ct) 


r=0 


3    cos    cot  /  _    3co   cos    cot    ,    sin   4cot    cos    cot 


+ 


)t    ( 
"9F 


TF 


sin      cot 
5    sin   cot    cos    2cot         sin   cot    cos    ^cot         cos    cot    sin    2cot 

2n 


ag(t) 


r=0 


6    cos    cot    /  3co    ,    sin    2ootN 

.    5     ^      W       3 ^ 

sm      cot  '  -^ 


a^(t) 


r=0 


3    cos    cot    /      sin   2cot         sin    ^cot         3co   cos    cot 

TZ —  ^ K  1^  8" 


sin   cot 


2 


+ 


sin   2cot    cos    ^(ot         cos    2cot    sin   ^cot         sin   cot    sin    2oot 


■9F 


96 


12 


cos    cot    sin    2cot  N. 


agCt) 


r=0 


P    cos    cot    /Sln   cot    _    sin   cot    cos    2cot 
5  .^         ^12  25 


sin   cot 
sin   cot    cos    ^cot         sin   cot    sin   ^oot         3^   cos    cot 


~w 


'W 


By 


.    cos    cot    sin    2cot    ,    cos    cot    sin   ^cotv 
■"  3 "■  9^ ) 


to   get    the    a.(t)   with   r  7^   0      let    3   ->   6+    ^^  cos   tot 

for   1   =   1,2,3    5      and      3    cos   cot   ->   3    cos   cot   +   ^ 

2m 


for  1   =   5,6,7,8, 
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b-,  (t)    =   b^(t)    =    C-,  (t)    =   ^^—  sin   cot    -   77  t    cos    wt 

1  J  1  <i(jO  ^ 

11  1  1 

b^(t)    =   bi,(t)    =    c^Ct)    =    T-  t    -   ^  t    cos    (jot+   Ti —  sin   2cot-   ^^r— sin   cot 

^         4         d  d  d  ^OJ  c:CO 

C.   N-Body  Renormallzation  Procedure 

It  is  to  be  observed  that  expression  (6.30)  has  a 
structure  somewhat  akin  to  equation  (5 -14)  and  therefore 

it  is  reasonable  to  expect  that  a  similar  method  for  renormal- 

t 
ization  of  divergences  can  be  employed.    In  actuality, 

however,  renormallzation  will  not  be  accomplished  by  way  of 

the  frequency  shifted  harmonic  terms  but  instead  by  their 

square .   At  first ,  it  would  appear  that  such  a  procedure 

alters  drastically  the  variationally  derived  equations  of 

motion,  but  it  will  be  shown  that  the  only  change  is  a  scaling 

of  frequency.   The  integrand  of   the  inequality  is  of  the  form 

(F  -F  ) 
(F^-AF)  e   ^   ^ 


where 


Fp^   is  the  basic  harmonic  oscillator  term 
AF   is  any  perturbation  term 

F-,   is  a  shifted  harmonic  term. 

2 
Furthermore,  by  choosing  G   -  F-,  one  can  write 


Fq-G 
e 


2 


2   F 
2iXG   -X     0  ,, 
e      e    e    dA 


and  so  one  has. 


Our  initial  approach  was  to  integrate  Inequality  (6. 30), 
make  the  complex  time  extension,  let^time  go  to  infinity, 
and  then  balance  the  various  secular  terms  by  their  orders 
of  infinity;  this,  however,  was  too  tedious  a  procedure. 
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(6.36)  e 


Fq-G 


exp  - 


2m3 


ll 


t  t 


0  0 


5(k,t' )C(-k,t")  dt '  dt" 


-  ir 


t 
r 


0 


x(k,t)^(-k,t' )  dt'  -  J-  r^  x(k,t)x(-k,t)] 


k 


+  2TTlk/M   2 


t 
/ 
k  0 


(JO  x(k,0)x(-k,0)   i  ^  /  Uk,t' )x(-k,t' )dt 


t  t 


I  H-t-  '^2 


•  e 


-{I   I    f    (-f2  )C(k,t')x(-k,t")dt'dt"}' 
k  0  0 

t  t ' 
+  i  ^  /  /  w^(k,t' )x(-k,t")dt'dt" 
k  0  0 


t  t 


■2  ^2 


(6.37)  = 


2iX{^  /  /  (to^-f]'')5(k,t' )x(-k,t")dt'dt"} 


dX  e 


k  0  0 


-X' 


-mB);  e 
k 


+  27Tlk/M   2 


t 
/ 
k  0 


w  x(k,0)x(-k,0)   i  ^  /  5(k,t' )x(-k,t' )dt 


•exp  -  75—5-  y  [ 

'^        2mB  f 
t 


t  t' 

r 


-  ir 


0  0 


C(k,t' )C(-k,t")  dt'  dt" 


1   2 


x(k,t)C(-k,t' )  dt'  -  ^r  x(k,t)x(-k,t)] 


2       ~2 
where   Ico   =  21Xa)  .   After  carrying  through  the  variational 

technique  one  obtains  the  frequency  scaled  equations 

I    x(k,T)  +  2X(w^-n^)x(k,T)  =  0 
I    ^'(k,T)  +  2X(aJ^-f^^)x(k,T)  =  0  , 


and  clearly  the  only  change  is  a  scaling  of  the  frequency 
which  leaves  the  calculation  essentially  unaltered.   Following 
the  methodology  of  Chapter  V  Section  D  one  can  again  show 
that  the  frequencey  shift  term  is  zero   and  the  only  finite 


t 


See  Appendix  D, 
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contribution  to  the  integral  comes  from  various  fourth 
power  terms.   This  leads  to  a  relatively  simple  inequality; 
hence,  following  the  above  outline  one  obtains  upon 
integration   the  finite  contribution  (for  expression  (6.30)) 


3tt 


(6.38)  li-   ^)6'|^(s,j,k,k,k,k)[ 


[c^-^(k)c22(k)-c^2^k)]^/2   ^^ 


1 


c-|^2*^^) 


c^^(k)    [c^^(k)c22(k)-cj2^1^^^^^^  ^"^  ^11^^^ 


c-|^2(^) 


3 


1 


[c-^-j^(k)c22(k)-c-j^2^1^)  ^    ^  c-j^-j^(k) 


2c-^2^^) 


c^^(k)[c^^(k)c22(l^)-Ci2(l^)^^''^  ^"^    ^^ 


2    2 
9   .    .,  X   roo  cos  cot 


2sin  cot 


^cl^ik. 


c^^(k)[c^^(k)c22(k)-c^2^^)^^''^  ^''   ^^ 


2    2 
9  .    .,  N   rco  cos  cot  ^  -, 
{c-,^{k)-   T )J 


2sin  wt 


r=0 
t->iT 


(6  3Q)  =  y     3tt3'['(s,j  ;k,k,k,k)a)  cosh  cot 


k  ^„3  3  6  .  ,  2  ,      2TTik/M^  cosh"coTx3 
23  m  CO  smh  coT(-e      +         ^ 


smh   ojT 


,2    2    4 


r/B"m"co'    r      2Trik/M,    cosh   cox         cosh'  cot -.x 
l{- ^ — [-e  +   T o — J) 


■1/2 


slnh   COT 


2  2 

sinh   ooT         sinh   cox 


,2    2    4    cosh    COT  /  3   m 


2    2    4 


m   CO 


c- 


4 — '  2 

sinh   OJT    sinli    wt 


oj      r_    27Tik/M      cosh   cot    _    cosh   cot-,-. 


-3/2 


2  2 

sinh   COT        sinh    ojt 


4           2    2    4                                                                 2 
,„4    4    8    cosh   COT, 3   m   oo      r      27Tik/M,    cosh   cot        cosh   wt-,  n -5/2-, 
+  3   m   to      n — ( T — L-e  +   o —  - — -i )  J 


smh   COT    smh   cot 


2  L 

sinh    (OX         sinh^'coT 


Y      3TT3t|^(s  ,  j  ,k,k,k,k)    cosh   cot 

k    ,,„4    4    6    .    ^4       ,      27rik/M  _^   cosh^coT.2 
43   m   CO   smh   ooT(-e  +   2 — ^ 

sinh   OJT 


t 


See    Appendix   D 
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[(- 


1 


sinh   wT 


r       2TTlk/M         cosh    OJT         cqshjiVT-,  s -3/2 

L  — e  +    2         ~  2 

sinh   GOT        slnh   cox 


.    2    cosh   (jot    /         1 

+    5 ( 2 

slnh   ojT        slnh   ojx 


r_    2TTik/M         cosh    cot         cosh    a)T-[s-5/2-| 
sinh   (jjT        slnh   ojt 


where    the    c.  .(k)    are    defined   in   Appendix   C    and   as    t   -^   °° 
this    leaves 


(f,   hn)    Y    3ttBi|^(s,j  ,k,k,k,k)      . 
^^•^°^    f    ^A   4   6.'      2wlk/M,3^ 


k    26-^m-^.-(l-e--^/^3^(l+e2-i^/")l/2 


(,^^2.1k/M)3/2 


1 n      Y    3tt3tj^(s   j,k,k3k,k) 

,.^    2^ik/M.5/2^~    f    ,,0^1    4    b",       27Tlk/Ms2 
(1  +  e  )  k    43   m   00    (1-e  ) 


[- 


+ 


(^^g27ylk/M)3/2  ,,^^^2wlk/M)5/2 


since    11m   tanh   cot   ->■    1    and    11m   1/sinh   cox   ->-   0 . 


Thus,    the    final    form   is 


t 


x-^o 


(6.41) 


dr 


X    e 


-B 


w^0+    -   3r 
r=0 

X->-oo 


A(B-B+1)    e 


-B 


r=0 

X->oo 


^    V    3Tnl^(s,j  ,k,k,k,k)  r 

f    ^    4,3   6/.      2Tflk7M73 


k    2m-33.6a-e2-ik/M^3    ^  (,,^2.ik/M)l/2         (,,32Trik/M)  3/2 


■  1 1   +   y    3^4^(s,j  ,k,k,k,k) 

,.^    2TTlk/M.5/2^        f     ,,    4_3    6    ,       2-^ik/Ivis2 
(1+e  )  k    4m   3    w    (1-e  ) 


^(l+e2^i^/M)3/2    -    (,^^2.ik/M)5/2^ 


This  is  so  because  X  was  taken  as  e''"—  exp-  p-  J^    r[x .  (t  )-x  (0 )  ]' 
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rfi   a?^    -   Y    3Tr4^(s,j  ,k,k,k,k) 


k   128m6^Yi    sin^^k/M      (l-e^^^^/^)^ 


(^_^2ulk/M)  (^^^27rik/M)l/2 


1 

2TTik/M. 


(1-e— -)         (1+3^-^^/^3/2 


(^_^2wik/M)  (i^^2TTik/M)5/2 


2(^^^2TTik/M)3/2         (^^^2ulk/M)5/2- 


where   y-i    =   Tiax    (A    .  ) 
1  sj 
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CHAPTER  VII 

Throughout  the  development  of  the  paper  the  scheme 
has  been  to  determine  the  transport  coefficient  given  the 
correlation  function.   This  was  accomplished  essentially 
by  the  method  of  Kubo  transforms.   We  now  consider  the 
hydrodynamic  point  of  view.   In  order  to  do  this  we  first 
review  a  paper  of  R.  D,  Mountain  [20]  where  he  calculates 
the  spectral  distribution  of  light  scattered  by  density 
fluctuations  in  a  dense,  monatomic,  one-component  fluid 
from  the  time  dependence  of  the  density  fluctuations 
predicted  by  the  linearized  hydrodynamic  equations  of 
irreversible  thermodynamics. 

A.   Hydrodynamic  Theory  of  Fluctuations 

Information  about  the  density  fluctuations  is 
contained  in  the  generalized  structure  factor  S(k,a))  which 
is  the  space  and  time  Fourier  transform  of  the  two-body 
correlation  function  of  the  medium.   The  decay  in  time  of 
density  fluctuations  is  presumably  described  by  the 
linearized  hydrodynamic  equations  of  irreversible  thermo- 
dynamics, which  are  the  continuity  equation 

Sp-, 

(7.1)  j^   +  Pq  div  V  =  0 

the  Navier-Stokes  equation 

2  2„ 

(7.2)  p      dv/dt   +   —  grad   p-,+   —  grad  T   -(^  n   +rig)grad  div  v=0 
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and  the  energy-transport  equation 

9T     c  (y-1)  3p.      p 
(7-3)      Po%  9-r  -  -^S 8t^  -  ^^  ^1  =  ° 

In  these  equations,  p  =  Pr,+P-|  is  the  number  density, 
T  =  T„+T-,  is  the  temprature,  p„  and  T^.  being  the  equilibrium 
values.   The  ratio  of  the  specific  heat  at  constant  pressure 
to  the  specific  heat  at  constant  volume  is  given  by  y=c  /c  . 
The  shear  and  bulk  viscosities  are  n   and  ripj  X  is  the  thermal 
conductivity,  6  is  the  thermal  expansion  coefficient,  and  Cq 
is  the  low  frequency  limit  of  sound  velocity.   For  this 
theory  of  small  oscillations  about  equilibrium,  the  transverse 
part  of  the  velocity  is  not  coupled  to  the  density.  Thus,  one 
can  neglect  the  transverse  part  of  the  Navier-Stokes  equation 
in  all  considerations.   This  limits  the  usability  of  this 
theory  to  fluids  for  which  angular  correlations  between  the 
molecules  are  not  important. 

To  proceed  one  first  eliminates  the  velocity  by  taking 
the  divergence  of  each  term  in  equation  (7.2)  and  then 
substitutes  for  div  v  from  equation  (7.1).   Next  one  obtains 
the  Fourier  (space)  and  Laplace  (time)  transforms  of  these 
equations.   Then  one  solves  the  resulting  equations  for  n(k,s), 
the  Fourier-Laplace  transform  of  the  density.  The  inverse 
Laplace  transform  of  n(k,s)  is  the  time-dependent,  k-th 
Fourier  component  of  the  density,  n(k,t): 

CO 

r 


(7.4)     n(k,s)  = 


dr 
v    0 


,,   -Ikr   -st    /   ,  V 
dt  e     e    p-,(r,t) 


-79- 


(7.^)'   T(k,s)  = 


dr 


dt  e"^^^  e"^^  T^(r,t) 


V    0 
The  elimination  of  the  longitudinal  component  of  velocity 
and  the  transformation  of  the  density  and  temperature  leads 
to  two  simultaneous  linear  equations  for  n(k,s)  and  T(k,s) 
in  terms  of  the  initial  values 


(7.5) 
(7.5)' 


n(k)  = 


T(k) 


dr  e  ^^^  P-,^(r,0) 


V 


dr  e 


-ikr 


T^(r,0) 


These  equations  are 

(7.6)  n(k,s)[s^+  Cgk^/Y  +(|  ng+ng)k^s/pQ]+(cQBpQk^/Y)T(k,s) 


4        2 
=  n(k)[s+(2-  ng  +  n3)k  /Pq] 


and 


(7.6)'  n(k,s)[-sc^(Y-l)/3]+  T(k,s)[pQC^s  +  Xk''] 

=  -n(k)  c^(y-1)/3  +  PqC^  T(k)  . 

The  solution  of  equations  (7.6)  and  (7.6)'  for  n(k,s)  is 

s^+(a+b)k^s  +  abk^  +  c^(l-Y~^)k^ 

s  +(a+b)k  s  +(cQk  +abk  )s  +  acQk  y~ 


(7.7)  n(k,s)  =  n(k) 


where   a  =  >^/PqC^  and  b  =  (^^  n^  +  ng)/PQ. 

The  next  step  is  to  compute  the  inverse  Laplace 
transform  of  n(k,s).  This  inversion  requires  that  one  find 
the  roots  of  the  denominator  of  equation  (7.7),  so 

(7.8)   s^  +  (a+b)k^s^  +  (cQk^+abk^)s  +  ac^k^Y""^  =  0 
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2 


f  2        2 

the  solutions   of  which  to  lowest  order  In  ak   and  bk   are 

s  =  +  Cpk  -  I  (a+b-aY~^)k^ 

2 
s   =   -(a/Y)k   ; 

or  in  terms  of  the  original  parameters 

(7.8)  s  =  +  ic„k  -  i  [  3  "^s  '^B  +  1_  (A_  _  A_)^^ 

-    °     2       Pq      Pq   c^    Cp 

(7.8)'  s  =  -  Xk^/p^c   . 

The  real  part  of  the  first  two  solutions  is  effectively  the 

low-frequency  sound-absorption  coefficient  and  henceforth 

2 
is  denoted  by  fk  ;  viz., 

(7.9)  rk^   =  h^-Il^J^  +  J^   (^  .  ^-)2k^    , 

2    PQ      PQ   %   ^p 

Using  thie  first-order  solutions,  equations  (7.8) 
and  (7.8)',  the  lowest  order  terms  of  the  inverse  Laplace 
transform  of  n(k,s)  are 

c  — c  2      c 

(7.10)  n(k,t)  =  n(k){^ — ^  exp  [ ^^  t]+  -^  exp[-rk^t]cos  c„kt} 

^p         PQ^p     ^p  " 


t 
In  R.  D,  Mountain's  paper  he  describes  a  convergent  scheme 

to  approximate  the  solutions  to  the  dispersion  equation  (7.8) 

by  using  a  power  series  with  parameters  ak  /cgk  and  bk  /cgk. 

These  are  ol'  the  order  of  1/100  because  typically  in  light 

scattering  one  has  changes  in  the  wave  vector  k   of  about 

lO-^cm"   (this  corresponds  to  a  scattering  angle  of  60°  when 

o 

6328A  He-Ne  laser  is  used  as  the  light  source). 
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The  k-th  component  of  the  density-density  correlation 
function  F(k,t)  Is 

(7.11)   F(k,t)  =  <n(-k)n(k,t)> 

c  — c  2      c 

=  <n(-k)n(k)>{^E — ^  ^^^^ Ak_  ^-j_^  _v  gxpE-fk^t  ]cos  c„kt} 

^p        ^O^p      ^p  ^ 

The  angular  brackets  <...>  Indicate  an  ensemble  average 
over  the  initial  values  of  the  variables.   The  quantity 
which  is  of  Immediate  interest  in  a  light  scattering 
experiment  is  the  generalized  structure  factor  S(k,co); 


00 

f 


dt  exp  (iwt)  F(k,t) 


0 


(7.12)  S(k,cj)  =  2  Re 

It  follows  from  equation  (7.11)  that 

S(k,a))  =  s(k)  a(k,w)  , 
where  S(k)  =  <n(-k)n(k)> 

is  the  ordinary  structure  factor  and 

2 
c  -c    2Ak  /pr.c 

(7.13)  a(k,a))  =  (-^y^)  g   °  P  g   g 

p    (Ak  /PqC  )  +C0 

+  ^  ^  rk^     ,      rk^ 


^p   (rk^)^+(cj+CQk)^   (rk^)^+(a)-CQk)^ 

B.   Thermodynamic  Aspects  of  Light  Scattering 

The  significant  parameters  in  the  above  a(k,ui)    can 
also  be  determined  by  quite  direct  thermodynamic  arguments. 
Inasmuch  as  the  scattering  of  light  is  directly  due  to 
density  fluctuations  (the  index  of  refraction  being  a  function 
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of  the  density  alone)  and  tlie  variation  of  density  is  due 
both  to  a  variation  of  the  pressure  and  of  the  temperature 
(or  entropy),  the  light  scattered  in  a  given  direction  owing 
to  the  density  waves  must  consist  of  two  parts:   a  mechanical 
[22]  one  with  an  altered  frequency  and  a  thermal  part  which 
must  have  the  same  frequency  as  the  incident  light.  The 
basic  ideas  of  this  theory  of  light  scattering  go  back  to 
the  work  of  Landau  and  Placzek  who  showed  it  to  be  especially 
appropriate  to  the  case  of  liquid  bodies  where  the  entropy 
fluctuations  are  much  larger  than  in  solids.   The  frequency 
of  the  light  scattered  by  the  fluctuations  at  constant 
entropy  (i.e.  phonons )  is  shifted  by  an  amount  proportional 
to  tl:ie  velocity  of  the  phonons.   Energy  and  momentum  considera- 
tions require  that  the  proportionality  constant  be  the  magni- 
tude of  the  change  in  the  wave  vector  of  the  scattered  light. 
Two  lines  are  observed  because  scattering  can  occur  from 
waves  traveling  in  opposite  directions  but  at  the  same  speed. 
These  lines  ar'e  broadened  somewhat  because  of  the  dissipative 
processes  which  damp  out  the  waves.   These  lines  are  known 
as  the  Mandelshtam-Brillouin  doublet. 

The  ratio  of  the  intensity  of  this  unshifted  line  1„ 
to  that  of  the  two  shifted  components  1-,  must  be  equal  to  the 
ratio  between  the  square  of  the  thermal  part  of  the  density 
fluctuation  Ap  =  Op/3T)   AT  or  Ap  =  Op/3s)   As   and  the 
square  of  its  mechanical  part  Ap  =  (3p/3p)rp  Ap  or  Ap=(g^)   Ap . 
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Following  the  method  of  Landau  and  Placzek  the  pressure 
and  the  specific  entropy  are  chosen  as  Independent  variables 
(which  is  a  reasonable  choice  for  relatively  long  waves 
capable  of  effectively  scattering  light)  and  the  volume 
V  of  a  given  small  mass  of  the  body  as  the  dependent 
variable,  in  place  of  the  entropy.   Then  one  will  have 

where  K_  is  the  adiabatic  bulk  modulus  and 

the  latter  relation  following  from  the  differential 
expression  for  the  total  heat 

(7.16)  dJ   =   T  dS  +  V  dp  . 

Denoting  the  minimum  value  of  the  latter  for  the  entire 

system  (including  the  incremental  volume  portion)  by  J„ 

we  have  further 

2 

(7.17)  -(AJ)3  ^  I  ^  (Ap)^  =  -  |(|J)3(Ap)^=  I  VK3(Ap)^  , 

9p 
and 

(7.18)  .(AJ)   %  1  ^  US)'    -    |(§)p(AS)^  =  jf  (AS)^  . 

Therefore  it  follows  that 

/(Ap)2  e-^^T)s/l<T 

^^•19)  (Ap)  =   ^   -(AJ),/kT =  79V-  =  —  ^S 

/  e  d(Ap)         ^ap-* 

and 
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P   -(AJ)  /kT 
/(AS)^  e      P     d(AS) 


/(e     P    cl(4S)        '•as'p 

—  00 

Thus, 


(7.21)  Iq/I,  =  (U>^(AS)^  /  (f  >s'(Ap)^  , 

and  using  the  Identities 
,3V.  /,9V,        3p  9T.  /.3T.       ,3S. 


and  the  thermodynamic  relation  (3V/3S)   =  (3T/9p)o   the 

p  o 

ratio  can  be  rewritten 

(7  22)  I  /I   =  (^)     (^) 


Moreover, 


^W't         ^dv'S         ^3sV^3V'' 


so 


T  (  ^  ^ ")       (  ^P  "\  _  (' ^P  "I      K    -K 

(7  97)   I£  =  rriP^  .rlE^  1  ^p  ^  ^  ^3v^T  ^3v^s  _  s  ^t  _  ^p  ^v 

^'•^^  I,    L^3V^T  ^3V^S^  ,3Sx        ,3^s  K„       c 

where  Krp=  -v(3p/3V)m  is  the  isothermal  bulk  modulus  and 

c   and  c   are  the  specific  heats  at  constant  pressure  and 
p       V  ''^  "^ 

constant  volume . 

Landau  and  Placzek  have  observed  that  the  widths  of 
these  lines  is  determined  by  the  lifetimes  of  the  density 
fluctuations  which  can  be  described  by  the  linearized 
hydrodynamic  equations  of  irreversible  thermodynamics. 

There  have  been  numerous  attempts  to  observe  the 
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spectrum  of  the  light  scattered  by  density  fluctuations.  It 
has  been  possible  to  detect  the  presence  of  Brillouin  lines  in 
liquids  using  the  techniques  of  high  resolution  spectroscopy. 

It  can  be  shown,  by  adapting  Van  Move's  neutron- 
scattering  results  to  light  scattering,  that  the  intensity 

I'  of  light  scattered  from  a  fluid  is 

2  4 
(7.24)         I'(R,a))  =  "  "^^g  Iq  sin^4.  s(k,a3)  . 

2ttc  R 

The  light,  scattered  at  the  origin,  is  observed  at  R.  The 
angular  frequency  of  the  scattered  light  is  Q, ,    u)  is  the  shift 
in  the  angular  frequency  and  k  is  the  change  in  the  wave  vector 
of  the  scattered  light  from  that  of  the  incident  light  in  the 
medium  k  •  k=k, 2  sin(0/2),  where  0  is  the  scattering  angle. 
There  are  N  spherically  symmetric  molecules  of  polarizability 
a  in  the  scattering  volume.  The  incident  light  is  assumed  to 
be  in  the  form  of  a  plane  polarized  monochromatic  wave.  The 
angle  between  R  and  the  electric  vector  of  the  incident  wave 
is  denoted  by  cj) .  The  polarization  of  the  scattered  light 
depends  only  on  (p   because  the  molecules  are  isotropic. 

Thus  it  is  seen  possible  to  study  the  spectral  distri- 
bution of  light  scattered  by  density  fluctuations.  Moreover, 
a  way  is  now  available  to  measure  transport  properties  of  the 
fluid  without  the  imposition  of  a  gradient  as  for  instance  the 
bulk  viscosity  can  be  inferred  from  width   of  Brillouin  lines. 


■4* 

The  width  of  the  central  line  is  proportional  to  the  thermal 
diffusivity,  X/qqC^.    Interest  in  this  quantity  is  greatest 
near  the  critical  point  where  conventional  measurements 
of  the  thermal  conductivity  are  difficult. 
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Appendix  A:   Crystal  Lattice  of  Coupled  Oscillators 

The  basic  Hamiltonian  is 

(A.l)      H  =  i-   y   pf  +   y   i-  q.A.,  q,   . 

j=l   '^    i,j     '^  ^ 

The  A.,  matrix  depends  upon  the  force  constants.  By 

transforming  to  normal  coordinates  {P^^}  and  {Q^^  via 


(A.2)      p.  =  I    c.^P^  ;      qj  -  I  o.,Q, 


k 
one  obtains  a  new  Hamiltonian 

(A. 3)     H  =  1  j^  (1  P^^  .  ..2  q2) 

where  w.  is  the  j-th  normal  mode  frequency  and  the  matrix 

J 

c .,  satisfies 
3^ 

(^•^)  ]    ^Jk^j*£   =   ^k£    • 

The  variation  of  the  normal  mode  Q^ '  s  and  Pj^'s  with  time 


follows  from 


2 
m  Qy.   +   i^y.Qy.   =  0  ; 


as  expected  one  finds 

P,  (0) 
(A. 5)      Qj^(t)   =  -^ —  sin  to^t  +  Q^(0)  cos  a)^t 

k 

P  (t)   =   Pi,(0)  cos  oj^t  -  mo)^  Qi^(O)  sin  oj^t  . 
The  momentum  of  the  j-th  particle  at  time  t  is: 

Pj(t)  =  I  a.^^p^^(O)  +  I   bj^Qi^(O) 
where 
(A. 6)      a.^  =  I    c.^c*^  cos  o^^t;   b.^=  "Cj^^^^k  ^^'^  V 
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From  orthogonality  one  has 


y  a .  .  =  y  cos  to,  t 

dj    k 


Also 


(A. 7)  I    (aji^-Kmo)^)  ^b^^^)  =  1. 

In   an  n-dlmensional   cubic    lattice   the    c.,  's   are 

(A. 8)  c        ■-■   M""/^    exp    27Ti(    I      j„k^)/M 

"^  a=l 

where  j  and  k  are  vectors  with  components  (j'  }  and   ^^(^,^ 

each  k,j  runs  through  integers 

j^  =  1,2,3,. ..,M;     k^  =  1,2,3,. ..,M. 

If  n   =   1    (linear   lattice) 

/,    n^  M-1/2  o    .    -^'a^a  M-1/2  n    •    '^'l^l 

(A. 9)      c        =   M  exp    2tti   — ^     =   M  exp    27ti   — ^ 

i=l?3  M'  k=123  M 


(A. 10)    ^   M   ^/^    exp    2TTi   -^    -M   ^^^   exp    -2TTi    ^tt^  =    6 


and   orthogonality   becomes 

L 

Jl  1    1 

Therefore 

,       2TTij    k   /M  27rij    k   /M 

(A. 11)  pj^=^^i-  \'  'Ji^^^^'  \- 

In  addition  one  has  [21] 

r.    .^^  V      ^kJl      2TTin£/M      -2Trink/M  2 

(A. 12)  ^      ^p  e  '       e  =   ^^n    ' 

k,£      ^^  " 

Moreover   the   transformation   is   unitary   since 

«  „  1      27Tij^k^/M  ^      -2uiJ^k2/M    ^ 


(A. 13)    I    q,    q,    =  I  7^  e  -^    ^         Q^^^    ^  e  Qj^^ 


k,  ,kp         12        12        k-,         11 


and  finally  since  k.,  was  a  dummy  index, 

it  At 

I       (q   q    -  Q   Q   )  =  0  . 
j^    ^1  ^1     ^1  ^1 

In  addition  the  normal  mode  frequencies  of  an 

N  =  M  ,  n-dimensional  nearest  neighbor  lattice  are  to  be 

used  and  are 

(A. 14)   muj^  =  4  I      y      sin^  -^    ;         k.  =  1,2,3,.  ..,M, 

j=l  ^  ^ 

where  y  is  the  force  constant  in  various  directions  and 
m  is  the  mass  of  each  particle.   The  frequency  spectrum. 

for  a  linear  lattice  is 

-1/2 
(A. 15)  g(a))    =    I    i^l-i^^) 

mijdj       =      4y-, 

2  2    ^^1 

mto        =      4y      sin     — j^  as   N  ^  °° 

Here  g(a))  dw  is  the  fraction  of  frequencies  between 
0)  and  oo+dw  while  w^  is  the  largest  frequency. 
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Appendix  B:  Explicit  Integration  of  Initial  Velocity 
Dependence . 
Integrating  out  x  (0)  from   A (B-B+l)e~^, . .   can  be 
accomplished  before  transforming  to  normal  modes.  Here,  the 


relevant  part  of  the  kernel  Is 


•2 


(B.l) 


-  ^  r  Zx.(t)x,(0)   -  i  Emex.(O) 
dx,(0)  e   "^    j  J     J     e   "^  J    J 


-1  /  E  C.(t' )x. (O)dt' 


Expressed  In  vector  notation  [23]  (B.l)  Is 


-  i-  rx(t)-x(0)   -  ^  mBx(0)-x(0) 
(B.2)   I  dx(0)  e      '^      "  -  e      ^        " 


-1/  C(t')-x(0)dt 


Integrating, 

CB.3)    {^)'''"exp   -jijC 


t       t' 
f 


C(t'  )-5(t")    dt'    dt" 


0    0 


1/ 


1      2' 


-   ir   J    x(t)-|(t')    df    -   ^  r^x(t)-x(t)] 
0 


for  which  explicit  use  of  the  formula 

1 

■  i   V  C  r,Y  n 

-a    aB-^  B 

e 


\-\ 


T  E  V  c  oYo   1  Z  D  V 
e  e  IT  dy^ 


(/2¥)^(Det   C)-l/2  g   2     a      aB  3 


was  made 
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Appendix  C:   N-Body  Variational  Equations 

The  N-body  kernel  under  the  normal  mode  transforma- 
tion is 

(CD  exp  -  jiji  [| 


t   t 


k 


C(k,t')C(-k,t")  dt'  dt" 


0  0 


t 
f 


-ir 


x(k,t)C(-k,t' )dt'  -  ^  r  x(k,t)x(-k,t)] 


„r   +2TTik/M   2  ,,  ^v  .  ,  „s 
-m&l    e-  0)  x(k,0)x(-k,0) 

k 


i  ^  /  ?(k,t' )x(-k,t' )dt' 
k  0 


i  .  ^.  V   2 
m 


Iff   mw  C(k,t ' )x(-k,t")dt'  dt 


k  0  0 


,1/2 


subject  to  cj  =  2(Y/m)  ^   sin  irk/M. 

The  variation  of  the  kernel  is  made  by  taking 

6/6C(k,T)  and  6/6x(k,T)  with  t  =  t',t"  and  k  =  +  k  (since 

the  kernel  is  symmetric  in  k).   This  gives 

t 
(C.2)  ^  (-  -i  J  ^(k,t')  dt'  +  2^  x(k,t)  +  ix(k,T) 
~        0         T 

I  oj^x(k,t")  dt")  =  0 


+  i 


t   0 


I      (-i?(k,T)+  i 

k 


to 


■C(k,t' )dt' )  -  0 


T-dif ferentiating  (C.2)  yields 

(C.3)  I    (x(k,T)-  co^x(k,T)  =  0  ;  I    (^(k,T)-  oj^?(k,T))  =  0; 
k    ~         ~  k    ~ 

and  putting  (C.3)  back  into  (C.2)  gives 
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(C.3)'  ^(C(k,T)  + 


T  t 

k 
which  reduces  to 


a3^5(k,t'  )dt'  )  =  0 
C(k,t' )  dt' )  =  0 


(C . 4 )    1    5(k,t)  =  0  ,   a  boundary  condition, 
k 


In  like  manner, 
t 

(C  .5)  li-      ^ 


m3 


0 


'  ~   C(k,t')dt'  +  2^  x(k,t)+  ix(k,T) 

T 


-  1 


x(k,t")  dt")  =  0 


implies 


(C.6)    I    C(k,0)  =  ^  J-  iroj^x(k,t)  +  im6w^x(k,0) 
k        "  k 

which  is  a  coupling  condition.   The  differential  equation 

set  (C.3)  when  solved  subject  to  (C.4)  and  (C.5)  and  the 

tautology  x(k,0)  =  x(k,0),  x(k,t)  =  x(k,t)  gives,  as  expected, 

the  solution 

x(k,t)    sin   oax   +    x(k,0)    sin    co(t-T) 


(C.7)      x(k,T)    = 


(C.8)       Uk,T)    =    IBmoo^C 


sm   ooT 
x(k,  t  )-x  (k  ,0  )cos    oat 


•    2    , 

sm   cot 


cos    a)(t-T) 


-,  p      x(k,t)cos    oat    -    x(k,0) 


+    p-  Iroa    ( 


•    2    , 
sm    cat 


)    cos    oa(t-T) 
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If  equations  (C.7)  and  (C.8)  are  now  reintroduced 

into  expression  (C.l)  it  becomes  after  algebraic  manipulation 

2  2  t 

-mB^  e^^iy^  '^   ^    (^0)   ir/^m3  I    x(k,t)  /  C(k,t')dt' 

(C.9)   e    ^  e  ^        "  ^        " 

-    P^^I-  F  r-^x^Ck.t))   ^  i  5;(?(k,t)x(k,t)-5(k,0)x(k,0)) 
^mpj^q      -      "^        k        - 
'   e  •<.  e    ~ 

which  after  further  arrangement  is  expressible  as 

exp   -   ^  I    (c^^(k)x^(k,0)+  2c^2(^)^(^'0)^^^'^)+  C22  (l^)x^(k,t ) ) 

with 

2        2                  2 
rw\    -    o^,  2    ^2Trlk/M        3mu   cos    wt    .    roo      cos    wt 
c-,-,  (.k;    =   pmw     e  +  p +   5 

sin   a3t  sin   cot 

2  2 

(ir^    =      3moo  rco      cos    oot 

sin  ut  sin   wt 

2  2 

/,  N              Bmco      cos    wt  1        2/1+cos    ojt  x 

^12^^  =  -  ^^ 2 -  2  ^"  ( 2 ^ 

sin  U3t  sin  tot 


"''see  Chapter  II,  equations  (2  .  ll)-(2  .  l4  ) 
for  an  analogous  calculation. 
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Appendix  D. 

I.   Typical  Multiple  Integrals  of  N-Body  Calculation. 

In  carrying  out  the  evaluation  of  the  three  third 

Yt "       3 
order  integrals  containing  e'   5(t')x  (t")  the  following 

typical  forms  arise  (for  which  terms  linear  in  y    are 

immediately  dropped  since  as  y  "^  0  they  contribute  zero). 


(DI.l) 


0    t 

r 


e^^"    sln^cot"    cos    co(t-t')    dt '    dt" 


_0O        —00 


yt  '    cos    (jj(t-t '  )  -3  „ 

P p [y    sin    ojt  '    -    3w    sin   cot'    cos    wt ' 

y      +    9oj 


)CjJ 


+    — p — — p-    (y    sin   cot  '    -    to    cos    cot')]    dt ' 
y    +    CO 


1    rcos    cot        00   sin   oot        co   sin   cot      cos    cot         sin   cot-i 
3C0    ^-^  „    2  „,  2  y  -^         203       J 


8 


CO 


32oj 


t      0 


(DI.2) 


t 
r 


e^^      sin^cot"    cos    aj(t-t')    dt '    dt" 


0        -00 


-O0J~ 


(y^  +  9aj^)(y^  +  oj^) 


cos    co(t-t'  )    dt ' 


3oo' 


sin   cot 
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II.   Representative  Quadratic  Square  Forms 


In  the  explicit  evaluation  of  inequality  (6.30)  one 

encounters  such  integrals  as 

7   -  k2c,5(k)x(k,0)x(k,t)  +  c„^(k)x^(k,t) 
(DII.l)   1-^1^  22 


—  00   _oo 


-  i(c^. (k)x^(k,0)         o 
e      -^-^  x(k,t)x-^(k,0)  dx(k,0)  dx(k,t) 


OO      CO 

r 


-00   _oo 


.  c   (k) 

.       -J   C5„(k)[x(k,t)+  ^      f^'.    x(k,0)]' 
x(k,t)x3(k,0)  e   2   22  C22(i<) 

-  ^  C   (k)x^(k,0)   i  /^,^v  x^(k,0) 
•  e   ^   ^^  e^  ^22^^^^         dx(k,0)  dx(k,t) 

Let  V  =  x(k,t)  +  ;; ^r^jY  x(k,0)    then  the  integral  becomes 


OO     OO 


c^^Ti<y 


[  (  ^  CnoCk)      4 


-  2^22'^^''^ 


—  00  _oo 


1         c  p(k)   2 
-  ^[c.,(k)-  /^.  v]x^(k,0) 
^   2   11      C22(k)  ^^  dx(k.O) 


-3ttc^2^^)^22^^^  /  2(c^^(k)c22(k)-c22(k))^'^2 

"   "   -  ^(2c,p(k)x(k,0)x(k,t)+c,,(k)x2(k,0)) 
(DII.2)   [   [  e   2    12  11 


.3, 


_oo  _co 


-  2    ^22  '^l^)x  (k,t ) 


-3tt  c^2^^)  *^11^^' 


2(c^^(k)c22(k)-c^2^^)) 


x(k,0)x^(k,t)  dx(k,0)  dx(k,t) 

-ryn      t)y  the  symmetric  exchange 


of  c^^Ck)  for  c-,-,(k).    The  :.^ecurrent  denominator  explicitly  is 


'22 


11 


(c-^-^(k)c22(l<)-C22'^^^-' 


t-^iT 


.„       2    2TTik/M       .       2    cosh    wx 
=    (  pmw   e  -    Bmiii      p — 

sinh   ojT 
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2    cosh    OJT  N  / 

-    roj      2 ^^- 

sinh   (jot 


3maj 


slnh   OJT 


2    cosh   OJT  s 
-    rco      2 ) 

sinh   OJT 


2  2 

/    o      2    cosh   OJT        roj(l+cosh   oot)n 

-    C-Bmoj      2 ^ 2 ^^ 

sinh   OJT        2    sinh   ojt 


and 


(c-^^(k)c22(k)    -    c-|^2(^)) 


2    2    4 
S   m   oj      ,      27Tik/M 


t^iT        sinh   OJT 
r=0 


cosh     OJT  cosh     OJT 

2        ~  2 

slnh   COT        sinh   ojt 


) 


Also  one  needs 


g^   (c-|^^(k)c22(l<:) 


2    /,  s  N-n/2 


r  =  0 

t^lT 


for  n   =   1,3,5, 


-n/2 


2 
(c^-|_(k)c22(k)-c-|^2^^^  ) 


(n+2)/2      |^(Cii(k)c22(k)-c^2^^)) 


r=0 

t^lT 


4  2 

n    gmoj   cosh   cjT/      2Trik/M      cosh   a3T         cosh   ojt  s. 

2       .    ,2  ^"^  .    ,  2         "       .    ,  2       '' 


sinh   OJT 


sinh   COT        sinh   a3T 


q2    2    4  „    .,    .^^  ,  ,2  (n+2)/2 

r  p    m   00         /_    2TTik/Ivi      cosh    OJT    _    c_osh_a)T_v -| 


smh   COT 


2  2 

sinh   COT         sinh   cot 
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III.   N-Body  Renormalization  by  Way  of  Square  of  Frequency  Shift 

As  mentioned  in  Chapter  VI,  the  renormalization 
procedure  for  expression  (6.30)  is  carried  out  about  the 
square  of  the  frequency  shift.   Typically,  one  has 


(DlII.l)     .  .  .m3(x(k,t)x(-k,t)-x(k,0)x(-k,t)cos    cot) 

+  27Tik/M 


ip(s,j,k,.  .  .)(b^(t)  +  l)+.  ..    =    [n^(e 


-1) 


•  (m3(x(k,t)x(-k,t)-x(k,0)x(-k,t)cos    cot)c    (t) 

+   m3(x(k,t)x(-k,0)    -   x(k,0)x(-k,0)cos    a)t)c2(t))]^ 

After  rearrangement, 

h 
(Dili. 2)  fi  =  {  .  .  .mB(x(k,t)x(-k,t)-x(k,0)x(-k,t)cos  tot) 

.,,/„     .    ,,            wSinojt         t    cos    cot     ,  \     /    f  r  r    —  ^^ 

•Wis  ,j  ,k,  .  .  .  )  {       2^ —   -    2 +...)/{L(e  -1) 

•  (m3(x(k,t  )x(-k,t  )-x(k,0)x(-k,t)cos    cot)(-   ^  t    cos    cot) 

+  m6(x(k,t)x(-k,0)-x(k,0)x(-k,0)cos  (jot)(|  t-  |  t  cos  cot)]^ 
+  linear  terms  +  oscillating  terms}. 

Now  the  numerator  either  oscillates  or  linearly  goes 

to  infinity  while  the  denominator   quadratically  goes  to 

2 
infinity  and  thus  fi   ^  0  as  t  ->•  °°.   A  second  characteristic 

form  to  be  examined  is 


(Dili.  3)  .  .  .[x(k,0)x(-k'  ,t)x(-k",t)x(-k'''  ,t)] 


2 


[a  (t)-  IBo^coso^t^ 
8y  sin  cot 


+  . 


[f22(e±2^i^/M_i) 


r=0 


(m3(x(k,t  )x(-k,t  )-x(k,0)x(-k,t)cos  ojt)c-,(t)  + 
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+  m3(x(k,t)x(-k,0)-x(k,0)x(-k,0)cos  wt )  C2(t)]^ 

where  all  the  a.(t)  terms  containing  a  l/y  contribution  have 
correspondingly  been  cancelled  in  pairs  (e.g.  a.,-*-  ag ,  ap-*"  &(-, 

ah^  ay,  a^->-  an)  before  letting  y  ^  0  because  all  are 

+ 
structured  alike   viz.  each  has  a  coefficient  like 

x(k,0)x(-k,t  )x(-k"  ,t  )x(-k '•' ,t )  which  is  just  a  fourth  power 

5  5 

of  position  times  a  3u)  cos  oot/sln  wt  or  gto/sin  cot. 

Hence  after  algebra  and  direct  substitution  for  c,(t), 

CpCt )  and  at- (t )  , 

(Dlll.ii)  ^^   =    .  .  .[x(k,0)x(-k' ,t)x(-k",t)x(-k"',t)] 

r„                ^^sin    4ojt    cos    wt    ,    5    sin   cot    cos    2cot 
•LB    cos    cot( n^ +   2T\ 

sin   cot    cos    4a)t         cos    wt    sin    2cot>-i, 
_^ ^ )J+... 

/   sin   ojt[m3x(k,t  )x(-k,t  )-x  (k,0  )x(-k  ,t )    cos    cot) 

•(2^  sin   cot-^-  t    cos    ojt)+   m3  (x  (k  ,t  )x  (-k  ,0  )-x  (k  ,0  ) 

•x(-k,0)    cos    cot)(77—  sin    2a3t    -   tt—  sin   a)t+   tt  t    -   77  t    cos    cot)] 
'  4a)  2w  2  2 

It  is  convenient  to  make  the  complex  extension  at 

this  point  and  so  one  may  consider  as  typical 

P  cr  -1 

(Dili. 5)       ...    sinh    4coT    cosh    cot/.  .  .  slnh    cot  (—  sinh    cot-t    cosh   cot) 

which    in   the    infinite    time    limit    goes    to    zero   and   once 

•       n2    . 
again   u      is    zero. 


f 

If  one  lets   t  ->  It   and  then   t  ^  «>  their  orders  balance. 
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Appendix  E:   Path  Integral  Linear  Terms. 

For  linear  terms  in  the  functional  Integral  one  has 
the  expression 


i  J 


-(y^Ay+B-y) 
ye       ~  ~   dy 


To  Integrate  such  terms  one  solves 


2  Ay  +  B  =  0 


•^2        ■^  ^stationary 

Let  y'  =  y  +  p-  A~  B  and  therefore 


1    1    -(y'^Ay')  -  \ 
(y'-  i  A  ^B)e         e 


-1, 


dy' 


T        1 
-y'  Ay'   -  ^ 

e        e 


'A  -'-B 


dy' 


since   y'  is  an  odd  function. 
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Appendix  F:   The  Reduced  Quartlc  Equation. 
The  quartlc  equation  [27] 

y   +  py   +  hy   +  ry  +  s  =  0 

may   be   reduced  to   the    form 

4  2 

z      +az      +6z+Y=0 

by  the  substitution  y  =  (z-  P/4 )  .   Let  Jl,m,n   denote  the 
roots  of  the  resolvent  cubic  t^+{^   a)t  +  (^^7   ^)t-  ^  B^=  0 
The  roots  of  the  reduced  quartlc  are 

z=  +(-/£-  /in-  /n)  ;    Zp=+(-/£+/m   +  /n) 
(F.l) 

z^  =  +  (  /I  -  v/^  +  /n)  ;    z^  =  +  (  /I  +  /m   -  /n) 

where  the  upper  sign  Is  used  If  B  Is  positive  and  the  lower 
sign  If  B  is  negative.  Therefore,  If  we  let 

c  =  z  -  p-  q     with   p  =  2q  ;    h  =  -  1 

s  =  q  ,        and    r  =  -  2  q  , 
we  have 

(z-  i  q)^+  2q(z-  |  q)^-  (z-  i  q)^-  2q(z-  ^  q)  +  q  =  0 

=  z^-  (1+  I  q2)z2+  (q3_q)z  +  (-  ^  q^  |  q^+  q)  =  0 

with  a  =  -(1+  I  q^);   6  =  q^-q;    Y  =  -  i|  q'^+  I  q^+  q. 

3     2 
The  cubic  equation   t  +a-,  t  fapt  +  a^  =  0  may  be  reduced  by 

the  substitution  t  =  (w-  a-, /3 )  to  the  normal  form  w  +6w+e  =  0 

12         13 
for  which  6=  ^rC  3a2-a-|^ )  ;  9  =  2T  (  3a^-9a.  a2+27ao  )  with 

solutions  £,m,n.  (i    =    A+B;   m,n  =  -(^+^')    +  i^  (A-B); 


A,B  =  [-  i  e  +  /e^/4  +  6^/27   ]^'^^ 


2 
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1  ct   -4y  6 

For   our    case   we   have    a-,  =    p-  a ;    a.^=   -^rr — -   ;    a^=   -   r~jr  and 

A  R  -  r    lr3c^^  +  9aB^    27  .2.    .1.  l,3a^  ,  9aB^ 


27  o2^t2 


2       3 

If  B  /4  +  6  /27  >  0,  then  there  are  one  real  root  and  two 

conjugate  imaginary  roots  for  the  cubic,  where 

6^  +  ^^  -     1    r  9a^  ,  Slg^g^  ,  (21.2.k     .  27     „4  2 

^  +  27  "  i^(27)2   ^^   (128)2  +  ^S:?^  ^  +  (8)(64)  «  ^ 

3(27)   ^3r2   9(27)  ^r4.   Ir   g^   ,  3a^Y    ,   Ba^y^   ,  y^n 
^"^^^^^       (64)^       2'  (48)3   4(48)2   (l6)(48)    ^^ 

Since  0  <  q  <  1  it  follows  that 
a^   =   -  (1  +  I  q2)3  <  0 

a3^  =   -  (1+  I  q2)q2(q-l)2(q+i)2  <  q 

B^   =   q2(q-l)2(q+i)2  >  q 

a^   =   (1+1  q2)2  >  0 

^4    3   2 
Y    =-i|-q+^q+q>0. 

By  examination  of  the  discriminant  near  the  end  points  of  q 
(and  at  any  intermediate  value)  it  is  indeed  found  to  be 
positive  definite.   Thus,  for 

4  % 
q%Y;^l;    3   '^0;     a;^^-2.5; 

2       3 
one  finds  the  leading  positive  term  of  9  /4  +  6/27  to  be 
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^ 

27  ^  (12)(64)^    27  ^100 


1  r-Jl^.Sr'  X  _  1  .  61 


1   5 

whereas  all  the  combined  negative  terms  are   -  py^Tno'^  ^^"^ 
therefore  the  discriminant  is  postiive  definite.   Likewise, 

for 

a,  J4  a,   4 

q'VY<<l;     3'^'q;     a;>^-l 

one  finds 

1  r  1       1        q   +     -1 

27  ^tfb"  "  110,592      3072  ^  •••  J 

which  is  also  positive  definite.  This  implies  that  the 

roots  of  the  quartic  are  complex  (i.e.  expression  (F.l)), 

and  hence  the  coefficient  of  the  t  term  cannot  reach  zero 
for  any  real  values  of  cos  ^q^- 
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